STABILITY ANALYSIS OF tt-KINKS IN A O-tt JOSEPHSON 

JUNCTION* 
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Abstract. We consider a spatially non-autonomous discrete sine-Gordon equation with constant 
forcing and its continuum limit(s) to model a O-tt Josephson junction with an applied bias current. 
The continuum limits correspond to the strong coupling limit of the discrete system. The non- 
autonomous character is due to the presence of a discontinuity point, namely a jump of tt in the 
sine-Gordon phase. The continuum models admits static solitary waves which are called 7r-kinks 
and are attached to the discontinuity point. For small forcing, there are three types of Tr-kinks. 
We show that one of the kinks is stable and the others are unstable. There is a critical value of 
the forcing beyond all static 7r-kinks fail to exist. Up to this value, the (in)stability of the ir-kinks 
can be established analytically in the strong coupling limits. Applying a forcing above the critical 
value causes the nucleation of 27r-kinks and -antikinks. Besides a 7r-kink, the unforced system also 
admits a static 3iT-kink. This state is unstable in the continuum models. By combining analytical 
and numerical methods in the discrete model, it is shown that the stable 7r-kink remains stable, and 
that the unstable jr-kinks cannot be stabilized by decreasing the coupling. The Sw-kink does become 
stable in the discrete model when the coupling is sufficiently weak. 
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1. Introduction. One important application of the sine-Gordon equation is to 
describe the propagation of magnetic flux (fluxons) in long Josephson junctions [18, 5]. 
The flux quanta or fluxons are described by the kinks of the sine-Gordon equation. 
When many small Josephson junctions are connected through the inductance of the 
superconductors, they form a discrete Josephson transmission line. The propagation 
of a fluxon is then described by the discrete sine-Gordon equation. For some materials, 
Josephson junctions are more easily fabricated in the form of a lattice than as a long 
continuous Josephson junction. In the strong coupling limit, a discrete Josephson 
junction lattice becomes a long Josephson junction. 

It was proposed in the late 70 's by Bulaevskii that a phase-shift of tt may occur in 
the sine-Gordon equation due to magnetic impurities [7]. Only recently this predic- 
tion is confirmed experimentally [35]. Present technological advances can also impose 
a TT-phase-shift in a long Josephson junction using, e.g., superconductors with uncon- 
ventional pairing symmetry [34] , Superconductor- Ferromagnet-Superconductor (SFS) 
TT-junctions [29], or Superconductor-Normal metal-Superconductor (SNS) junctions in 
which the charge-carrier population in the conduction channels is controlled [4] . 

A junction containing a region with a phase jump of tt is then called a O-tt Joseph- 
son junction and is described by a O-tt sine-Gordon equation. The place where the 
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0-junction meets the Tr-junction is called a discontinuity point. A O-tt Joscphson junc- 
tion admits a half magnetic flux (scmifluxon) , sometimes called n-fluxon, attac;lied 
to the discontinuity point [14]. A semifluxon is represented by a vr-kink in the O-tt 
sine-Gordon equation [32] . 

Using the technology described in [14], a O-tt array of Joscphson junctions can be 
created as well. Such system can be modeled by a discrete O-tt sine-Gordon equation. 
A short numerical study of a discrete 7r-kink is given in [31]. 

The presence of the scmifluxon in a O-tt .Joscphson junction or a O-tt array of 
Joscphson junctions opens a new field where many questions, that have been discussed 
in detail for the 27r-kink (fluxon) in the sine-Gordon equation, can be addressed for 
the TT-kink too. The fact that the 7r-kink cannot move in space, even in the continuum 
case, will give a different qualitative behavior such as the disappearance of the zero 
eigenvalue (Goldstone mode), as will be shown later. 

In this paper we will study both the continuous and discrete O-tt sine-Gordon 
equation, especially the stability of the kinks admitted by the equations. Knowing 
the eigenvalues of a kink is of interest for experimentalists, since the corresponding 
eigenfunctions (localized modes) can play an important role in the behavior of the 
kink [27]. 

The present work is organized as follows: in section 2 we will describe the math- 
ematical model of the problem and its interpretation as a Joscphson junction system. 
We will discuss the discrete system as well as several continuum approximations. In 
section 3 we consider the continuous O-tt sine-Gordon equation which describes a con- 
tinuous long Joscphson junction with discontinuity point. It is also the lowest order 
continuum approximation for the discrete system, not reflecting any lattice spacing 
(coupling) effects. In [32] it is shown that there exist three types of 7r-kinks in the O-tt 
sine-Gordon equation. We will analyze their stability and show that one type is sta- 
ble and the other two are unstable. A higher order continuum approximation, which 
includes terms representing a small lattice spacing (strong coupling), is considered in 
section 4 . It is shown that for small values of the lattice spacing parameter, the three 
types of TT-kinks persist and their stability properties do not change. In section 5 the 
discrete O-tt sine-Gordon with large lattice spacing (small coupling) is analyzed, es- 
pecially the existence and stability of 7r-kinks. Numerical calculations connecting the 
regions of small and large lattice spacing (weak and strong coupling) will be presented 
in section 6. In this section the analytical results of the previous sections are linked 
together. Conclusions and plans for future research are presented in section 7. 

2. Mathematical models for O-tt junctions. 

2.1. The discrete O-tt sine-Gordon equation. The Lagrangian describing 

the phase of a O-tt array of Joscphson junctions is given by 

2 1 / i i \ 2 



dt. 



where </>„ is the Joscphson phase of the nth junction. The phase jump of tt in the 
Joscphson phase is described by 9n where 



(2.2) 



0, n < 0, 
-TT, < n. 



The Lagrangian (2.1) is given in nondimcnsionalizcd form. The lattice spacing 
parameter a is normalized to the Joscphson length Aj, the time t is normalized to the 
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inverse plasma frequency oj^^ and the applied bias current density 7 > is scaled to 
the critical current density J^. 

The equation of the phase motion generated by the Lagrangian (2.1) is the discrete 
O-TT sine-Gordon equation 

(2.3) ^2 = ~ sm(0„ + 6i„) + 7. 

We use n e Z for the analytical calculations, but of course, the fabrication of the 

junction as well as the numerics arc limited to a finite number of sites, say IN . We 
will take the boundary conditions to represent the way in which the applied magnetic 
field h = H/{XjJc) enters the system, i.e., 

/2 (^-Af+i - (f)-N ^ (l)N - (t)N-l ^ ^ 

a a 

In the sequel we will always consider the case when there is no applied magnetic field, 
i.e., we will take h = 0. 

2.2. Approximations to the lattice spacing in the continuum limit. 

There are various continuum model approximations for (2.3) that can be derived 
in the continuum limit a <^ 1. Writing (?i„ = (f){na) and expanding the difference 
terms using a Taylor expansion gives 



and 



1 1 ^ h 

Z>7i+1 -<Pn \- a 



h—Ci ' ' 



fe=0 

Thus the continuum approximation for (2.3) is 

(2.5) ^tt - La(l>xx = - sin(0 + 6*) + 7, 

where ^(a;) is defined similar to (2.2), i.e.. 



e{x) 



0, a; < 0, 
— TT, a; > 0. 



The continuum approximation for the Lagrangian is 



dx dt 



Note that the normalizations in the discrete system imply that the spatial coordinate x 
is normalized to the Josephson length Aj. 

There are several ways to derive approximations for the operator La when a ^ 0, 
see for example [28] . The first obvious approximation is 

(2.6) (l)tt- (l)xx- r^4'xxxx = -sin{(j) + e)+'^, a; ^ 0. 
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Another approximation can be found by using that (1 — j2^x^) = ^ — ^d^x + — 
This results reflects the invertibility of La up to fourth order. Hence (1 — j2^^^) 
acting on (2.5) gives the approximation (up to fourth order terms) 

(2.7) <^xa= = </'tt+sin(0 + 6i)-7-— 5^^((/)u+sin(0 + 0)), a; ^ 0. 

Expanding this equation and using the expression for (j)xx again, we get 

(2.8) (t)xx = (t>tt+sm{(j) + e)-'Y 



2 

tttt 



12 

+ cos((^ + 6l)[?!>„+sin(^!> + 6')-7]), a; ^ 0. 



The steady state equation for (2.6) is 



hx + Y^(l>xxxx = sin((/) + 6*) - 7, Xt^O, 



while (2.7) yields the equation 



(1 - —dxx) sin((/) + 6*) - 7, a; ^ 0, 



and (2.8) gives 



a2 

^xx = sin(0 + 9)- J- ■^{-<t>l sin(</) + 9)+ cos(</. + 0)[sin(</. + 0) - 7]), x^O. 

Unfortunately the last two equations are not Hamiltonian, so we have lost the Hamil- 
tonian properties of the original system, while the first equation is singularly per- 
turbed. 

Yet another approximation that has a variational structure and is not singularly 
perturbed can be obtained by combining the two equations that have lost their vari- 
ational character. Indeed, taking (2.7) twice and subtracting (2.8) gives 

(2.9) 4)xx = 4>tt + sin{4> + 
a? 

- — {2(j)xxtt + 2(j)xx cos{(f> + 9)-(j)l sin(0 + 9)- (j)uu - (j^tt cos(0 + 9) 
+ ^1 sin(0 + 9)- cos{<p + 9){<ptt + sin(0 + 6») - 7)) , x ^ 0. 
The Lagrangian for this system is 

/ 5'^?-5'^^-l + cos(</' + ^)+7</' 

1 „1 

dxdt . 



a" 



,dx{(t)tt + sin(^ + e)) + ^{(l)tt + sin((/) + 9)- -ff 



The static equation for (2.9) is 
(2.10) = sin(0 + 61) - 7 



,2 



- — {2(t)xx cos((/> + 9)-(l)l sm{4, + 9)- cos(0 + 6i)(sin(0 + 6») - 7)) , a; ^ 0. 
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This equation is a regularly perturbed Hamiltonian system with the Hamiltonian 

2 2 

2(1 + ^ cos(^ + 61)) 24 
which implies p = (f)^ ( 1 + ° coh(0+e) 



6 

In this paper, we will analyze equation (2.9) as a continuum strong interaction 
limit which incorporates some effects of the lattice spacing into the model. The model 
equation (2.9) is chosen as it is non-singular and has the same conservative properties 

as the discrete system, reflecting its physical properties. 

3. The TT-kinks and their spectra in the continuum limit. In this section, 
we will consider (2.9) for a = 0, which is a model for an ideal long O-tt Josephson 
junction: 

(3.1) 4>u-(l>xx+s.m{(j) + e)=^, a; 7^0. 

For a Josephson junction without an applied bias current or a phase jump, i.e., for 
7 = and 0{x) = 0, the model corresponds to the sine-Gordon equation. A stable so- 
lution of the sine-Gordon equation is the basic (normalized) stationary, monotonically 
increasing fluxon, given by 

(3.2) (/>flux(a;) = 4arctane^, 4>^^^{Q) = -k 
(see [10]). 

In general the discontinuous function 6{x) in (3.1) will introduce a discontinuity at 
X = Q for the second derivative (i)xx- Hence, the natural solution space for (3.1) consists 
of functions which arc spatially continuous and have a continuous spatial derivative. 
The behavior at infinity is regulated by requiring that the spatial derivative of the 
solution belongs to Hi (M) (which allows the phase to converge to a nonzero constant 
at infinity). Therefore, the equation (3.1) is considered as a dynamical system on the 
function space 

M ={(/>: M ^ M I e i?i(M)}. 

It is straightforward to find that for I7I < 1 and a; < 0, the "fixed points" of (3.1) 
are = arcsin(7) and = — arcsin(7) -|- n. Similarly, for I7I < 1 and x > 0, 
they are (pf = arcsin(7) -I- tt and 0+ = — arcsin(7) + 27r. In [32], it is shown that 
there exists various types of stationary fronts, which connect the equilibria. Most 
stationary fronts are so-called 7r-kinks, which are static waves connecting equilibrium 
states at x = ±00 with a phase-difference of tt. Such waves are solutions of the static 
wave equation 

(3.3) (j>xx - sin((/) + 9) = -7, x ^ 0. 

In the x-dynamics of (3.3), the points (f>f are saddle points and the points with 0^ 
are center points. Thus a 7r-kink connects (f>J with (f)'^ . 

In this section we will consider the stability of those 7r-kinks. For completeness, 
we first describe the various types of 7r-kinks as found in [32]. These 7r-kinks are 
constructed by taking suitable combinations of the phase portraits for 6 = and 
9 = —TT. The phase portraits for 7 = are essentially different from the ones for 
< 7 < 1 (the case — 1 < 7 < follows from this one by taking cf) ^ —cj) and 
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Fig. 3.1. (a)The phase portrait of system (3.3) for 7 = 0. The trajectories for x < are 
indicated with bold lines, the trajectories for a; > with dashed lines. Any orbit of (3. 1 ) switches at 
.i: = from bold to dashed,. TJie type 1 semifiuxon switches at di and corresponds to one of the gray 
arrow-lines. The 3n-fluxon switches at d2 and is denoted by the other gray arrow-line, 
(b) The phase portrait of system (3.3) for ■y = 0.1. For simplicity, only the stable and anstable 
manifolds of the fixed points are shown. Apart from d\, there are also the points d2 or di which can 
be used for the switch position of x = to obtain a solution with a phase difference n between the 
end points. 



7 —7). In case 7 > there are homoclinic connections at kn + arcsin(7), k G Z, k 
even {9 = 0) or fc odd {9 = — tt). If 7 = 0, then these homoclinic connections break 
to heterochnic connections between kn and (fc + 2)?:. 

The phase portrait of (3.3) for 7 = is shown in Figure 3.1(a). Following the 
notation in [32], in case 7 = 0, there are two types of heteroclinic connections (kinks) 
in the O-tt junction. The first one, called type 1 and denoted by (j)\{x; 0), connects 
and TT. The point in the phase plane where the junction lies is denoted by c?i(0). The 
second one. called type 2 and denoted by (63^ (.x; 0), connects and 37r. Now the point 
in the phase plane where the junction lies is denoted by 1^2(0). This solution is not a 
semifluxon, but it will play a role in the analysis of some of the semifluxons for 7^0. 

If < 7 <C 1, then there are three types of vr-kinks (heteroclinic connections) in 
the junction, all connecting arcsin(7) and tt + arcsin(7). A phase portrait of (3.3) for 
nonzero 7 is shown in Figure 3.1(b). The first semifluxon, called type 1 and denoted 
by 0^(.7': 7), is a continuation of the connection at 7 = 0. The point in the phase plane 
where the junction lies is denoted by 1^1(7). The 7r-fluxon 0^(x;7) is monotonically 
increasing. 

The second one is called type 2 and denoted by 4>^{x] 7). In the limit for 7 ^ 0, it 
breaks in the 37r-kink and the heteroclinic connection between 2>n and tt (a — 27r-kink 
or an anti-fluxon) . The point in the phase plane where the junction lies is denoted by 
^2(7)- The TT-fluxon cfy^ix]^) is not monotonically increasing, but has a hump. 

The third one is called type 3 and denoted by 0^(x; 7). In the limit for 7 — » 0, it 
breaks in the heteroclinic connection between and 2^ (fluxon) and an anti-semifluxon 
like the type 1 wave but connecting 27r and tt. The point in the phase plane where the 
junction lies is denoted by ^3(7). This 7r-fluxon has a hump too, but a lower one than 
the type 2 wave. Following the first homoclinic orbit, the junction points are ordered 
such that (ii(7) comes first, followed by (^2(7), followed by (^3(7) (sec Figure 3.1(b)). 

If 7 increases, the points (^2(7) and (^3(7) approach each other, until they coincide 
at [32] 

(3.4) 7 = 7* = _L= 

in the point (tt + arcsin(7*), 0). At this point, the type 2 wave ^^(a;;7) ceases to 
exist (in the limit it breaks into half the homoclinic connection for x < and the 
full homoclinic connection for .x > 0). The type 3 kink (j)^{x\^*) consists of half the 
homoclinic connection for x < and the fixed point for a; > and this wave can be 
continued for 7 > 7*. For 7 > 7*, the type 3 kink is monotonic. 
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If 7 increases further, the points ^1(7) and ^3(7) approach each other [32] until 
they coincide at 

2 

(3.5) 7 = 7„ = -. 

TT 

When 7 = 7cr, the orbit homochnic to the hyperbohc fixed point for x < is tangential 
at di (7) = (7) to the non-homoclinic stable manifold of the hyperbolic fixed point 
for X > 0. As soon as 7 > 7cr, there is no more intersection of the homoclinic orbit 
for a; < with a stable manifold of the hyperbolic fixed point for x > 0. This implies 
that no static 7r-fluxons can exist. For more details, see [32]. 

After recalling the description of the 7r-kinks from [32], we can start the stability 
analysis. It will be shown that the type 1 7r-kink is nonlinearly stable for all < 7 < 
7cr. The type 2 and type 3 7r-kinks are linearly unstable for all values of 7 for which 
they exist. First wc consider the linearization about the 7r-kinks. 

Theorem 3.1. The linearizations about the various ir-kinks have the following 
properties: 

(i) The eigenvalues of the linearization about the m,onotonic type 1 n-kink cbl^ix; 7) 
are strictly negative for < J < jcr- At j = jcr, the largest eigenvalue is zero. These 
TT-kinks are linearly stable. 

(ii) The largest eigenvalue of the linearization about the type 2 n-kink (j)\(x;^) 
is strictly positive for < 7 < 7* . These -kinks are linearly unstable. 

(iii) The largest eigenvalue of the linearization about the type 3 it -kink (l)\{x;^) 
is strictly positive for < 7 < 7cr. These tt -kinks are linearly unstable. In the limit 
for 7^0 and 7 7cr, the largest eigenvalue converges to zero. 

REMARK 3.2. Note that the instability of the two non-monotonic 7r-kinks cannot 
be established by the classical Sturm-Liouvillc argument. In the classical, auton- 
omous setting, the derivative of the wave about which the system is linearized, is 
an eigenfunction of the linearized system. This eigenfunction is associated with the 
translation invariance of the original system and hence corresponds to an eigenvalue 
A = 0. If the wave is non-monotonic, then its derivative has a zero, which implies 
that A = is not the largest eigenvalue [33] and that the wave must be unstable. Due 
to the discontinuity at a; = 0, our system is non- autonomous, thus not invariant with 
respect to translations, and A = is (in general) not an eigenvalue. Thus, it cannot 
a priori be concluded that the non-monotonic 7r-kinks must be unstable. 

To prove Theorem 3.1, it will be shown that the linearization about a 7r-kink has 
an eigenvalue zero if and only if the 7r-kink takes a value which is a multiple of tt at 
x = Q. Since the value at a; = is related to the point ^^(7), it can be derived that this 
happens only at 7 = 7cr for the colliding type 1 and type 3 waves. To complete the 
proof, we will derive expressions for the largest eigenvalue of the linearization about 
each semikink near 7 = in three separate lemmas and use that the eigenvalues are 
continuous in 7 to derive the sign of the largest eigenvalue on the existence interval 
of the TT-kink. 

To linearize about a solution 4>\^{x\ 7), write (i){x, t) = (l)l^{x; 7) +v{x, t), substitute 
this in the model equation (3.1) and disregard all higher order terms: 

(3.6) [D^^ - cos{cl>i{x;^)+e{x))]v = Duv. 

Using the spectral Ansatz v{x, t) = e'^^v{x), where v{x) is a continuously diff'erentiable 
function and dropping the tildes, we get the eigenvalue problem 



(3.7) 



C{x;'y)v = X^v, 
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where D is defined as 



(3.8) r(x;7) =D,^-cos(<^;(a;;7)+%)). 

The natural domain for D is _ff2(lR)- We call A an eigenvalue of O if there is a 
function v S Hii^, which satisfies ^{x; = Av. Since depends smoothly on 7, 
the eigenvalues of will depend smoothly on 7 too. 

The operator is symmetric, hence all eigenvalues will be real. A straightforward 
calculation gives that the continuous spectrum of is in (—00, — y'l — 7^). Since 
the eigenfunctions are continuously differentiable functions in H2{M.) by the Sobolev 
Embedding Theorem, Sturm's Theorem [33] can be applied, leading to the fact that 
the eigenvalues are bounded from above. Furthermore, if vi is an eigenfunction of 
with eigenvalue Ai and V2 is an eigenfunction of with eigenvalue A2 with Ai > A2, 
then there is at least one zero of V2 between any pair of zeros of vi (including the 
zeros at ±00). Hence if the eigenfunction vi has fixed sign, then Ai is the largest 
eigenvalue of 

The following lemma gives a necessary and sufficient condition for to have an 
eigenvalue A = 0. 

Lemma 3.3. The eigenvalue problem 



has an eigenvalue A = if and only if one of the following two conditions holds 

(i) Dxx'P^ix]^) is continuous at x = 0, i.e., ^^(0;7) = kn, for some k G Z; 

(ii) -D2;0^(O;7) = and there are some x±, with sgn(a;±) = ±1, such that 
Dxcj>i{x±;-i)^Q. 

Proof Since (f)l^{x;j) converges to a saddle point for 00, this implies that 

£'s</'7r(a;; 7) decays exponentially fast to for \x\ — > 00. Since (^^(a;;7) solves (3.3), 
differentiating this ODE with respect to x, gives 

C{xn)DMxn) = o, for xj^o. 

This implies that for any constant K, the function w'j({x) = K Dx(l>\{x;^) satisfies 
0{x; 7) w^j^{x) = for a; 7^ 0. Hence for any K_ and Kj^, the solution 



w\x) 



w]^_ (x), X <0, 



solves £'(a;;7) w^{x) = for a; ^ 0. The function w^{x) is continuously differentiable 

if and only if the following two conditions hold 

1. w]^_iO-) = w],^{0+), in other words, D^(/.j,(0; 7) = K+ D^cj^HO;^), 
since (j)^^ is continuously differentiable; 

2. D,:w},_{0-) = a,w},^{0+), thus Z?,,<^;(0-;7) = K+ Z?,,,^;(0+; 7). 
The first condition is satisfied if A'_ = K+ or D.x<6;(0: 7) = 0. If Dx(f)l{{); 7) = 0, we 
can choose K± such that the second condition is satisfied and we do not end up with 
the trivial solution, except when D^cji^^ix; 7) is trivial for either a; > or x < 0. 

If Dj;(/)^(0; 7) 7^ 0, we need Dxx4>lr to be continuous at x = in order to satisfy the 
second condition. Since Dxx(f'lv{^^ l) — sin(0^(x; 7) + 0{x)) — 7, Dxx<Pv is continuous 
at X = if and only if sin((/)^(0; 7)) = 0. These arguments prove that if one of the 
two conditions are satisfied, then A = is an eigenvalue of £*. 



STABILITY ANALYSIS OF tt-KINKS 



9 



Next we assume that A = is an eigenvalue of hence there is some continuously 

differentiablc function such that £*(.t)?;*(x) = for x ^ and (/(x) for 

I a; I 00. The only solutions decaying to zero at +oo are the solutions on the one- 
dimensional stable manifold and similarly, the only solutions decaying to zero at — oo 
are the solutions on the one-dimcnsional unstable manifold. The stable and unstable 
manifold are formed by multiples of Dx4>\- So we can conclude that there exist K± 
such that 



K_Dx<p\{x) for a; < 0, 
K+Dx(l)i{x) for a; > 0. 



Now we are back in the same situation as above, so we can conclude that either one 

of the two conditions in the lemma must be satisfied. □ 

The second condition in the lemma does not occur. Indeed, the first part of the 
second condition, i.e., Dx<p)^{0;'y) = happens only if di has its second coordinate 
zero, hence only at 7 = 7* with ^2 = ^3. At this point, the solution 0^(2;; 7*) has 
ceased to exist and the solution (j)^{x; 7*) consists of the fixed point for a; > 0. Hence 
this solution does not satisfy the second part of the second condition. 

To see for which value of 7 the first condition is satisfied, wc derive the relation 
between (/'^^(O; 7) and 7. Multiplying the static equation (3.3) with D^cj)]^ and rewriting 
it gives 

i?,[(Z?,0;(.T;7))2] =2i?,[-7</);(a:;7)-cos(</.t(a:;7)+0(x))], x ^ 0. 
Integration from ±00 to and using that D2;(/)^(±oo; 7) — 0, shows 

iD,<PliO;j)f = 2[-7(</.;(0;7) - - cos(,^;(0; 7)) + cos((/.;(-oo; 7))], 

iD,^l{Orf)f = 2[-7(</.;(0;7) - ^U+oo; 7)) + cos(^;(0; 7)) - cos((/.;(+oo; 7))]. 

Subtracting these two equations and using that (j)l^{+oo; 7) = ^5r(~°o> l) + S^^ 
that 

(3.9) = -7r7-2cos(</.;(0;7)), hence cos(</>;(0;7)) = ^. 

Thus the first condition is only satisfied when cos((?!)^(0; 7)) = ±1, hence 7 = |- = 7cr- 
The following step in the analysis of the eigenvalues of the linearization is to 
consider the behavior of the eigenvalues for 7 small. First note that at 7 = 0, we have 
an explicit expression for the 7r-fluxon and the 37r-fluxon (see (3.2) for the expression 
of 0flux): 

('^^f)\ ^W^.m _ / <^flux(a;-ln(l + \/2)), for a; < 0, 

(3.10) cp^ix, 0) - I ^ _ ^ ^^^^ for X > 0, 

C^in ^2 _ / 0flux(a; + ln(l + ^2)), for x < 0, 

(3.11) <^3.la;,U)-| 

Hence the derivatives of both functions are even and cos(0^(x; 0) + 9) is continuous 
and even, since 0^(0; 0) = | and (/)§^(0; 0) = ^. 

For 7^1, the homoclinic orbit in the system with 9 = will be crucial for the 
approximation of type 2 and type 3 solutions. This orbit is homoclinic to arcsin(7) 
and will be denoted by cf)h{x;^). It can be approximated up to order 7 by using the 
27r-fluxon ^aux and its linearization. 
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Lemma 3.4. For 7 small, we have for the even homoclinic connection (j)h{x;'y) 
(3.12) Mx;i) = <l>fiuAx + L^il)) + 1 + Lr^il)) + l^R2{x + L^{j);-f), x < 0, 
where the expression for the 2TT-fluxon 0flux can he found in (3.2), 



— 1 + cosh X 



cosh^ 

and Lt^(^) is such that (/)/i(— ^^(7); 7) = tt 



1 



— arctan e 



/ X 



coshx 
= '?^'flux(0), implying 



cosh a; 



inha;^ 



+ sinh 



(3.13) 



L,(7) = i|ln7|+ln^ + 0(V7)- 



Furthermore, 7^i?2(a; + L7r(7); 7) = ^{l), uniform for x < anrf 7^1(1/^(7); 7) = 
0{^). Thus 



(3.14) 



f.^0) = 27r- 2^^^7 + 0(7) • 



Finally, 4>i{x;^) = 0(1) and R2{x;j) = 0(1), uniform for x < 0. 

Proof It is more convenient in the following perturbation analysis to follow the 
normalization of 4>Rux{x), i.e., in this proof we introduce new coordinates x = x+L,r(7) 
where L^{'y) is such that (^/i(— i7r(7);7) = tt = (^flux(O). In the following we will drop 
the tildes and work in those new coordinates. As (l)h in the original coordinates was 
even, we get in the new coordinates DxCphiLvil)',!) = 0. This condition will be used 
later to determine an asymptotic expression for i/7r(7). 

In the new coordinates, we introduce the expansion 

^h{x;'y) = 0flux(a;) + 70i(a;) + 7^i?2(a;;7), x < L^i-y). 

By linearizing about ^aux, it follows that the equation for (pi is 

(3.15) £.{x) (f>i = -1, where C{x) = D^x - cos((?!>flux(a;)). 

The operator C{x) is identical to the operator associated with the stability of (/'flux(a;). 
The homogeneous problem Cij} = has the following two independent solutions, 



(3.16) 



ipb{x) 



cosho 



1pu{x) 



cosh: 



- sinh a;, 



where '06 (a;) = ^^(/>flux(a^) is bounded and ipu{x) unbounded as x - 
variation-of-constants method, we find the general solution to (3.15), 



±00. By the 



hix;A,B) = 



4 1 u 1 
^ + - cosh X + - 

+ [B - arctan e^] 



cosh^ 



d^ 



1 



cosh a; 



cosh a; 



sinh X 



with A,B G R. The solution 0i(a;) of (3.15) must be bounded as a; ^ —00 and 
is normalized by ^!>i(0) = (since (phiO) = ^flux(O) = tt). Thus, we find that A = 
— i and _B = 0. Note that \imx^-ca 4'i{x) = 1, which agrees with the fact that 
lima;^_oo (ph{x) = arcsin7 = 7 + 0(7"^). The solution 0i(a;) is clearly not bounded as 
a; — > 00, the unbounded parts of <^i(a;) and ■^<pi{x) are given by 

(3.17) <l>i\u{x) = — arctan sinh a;, -^Muix) = — arctan cosh a;. 

dx 



STABILITY ANALYSIS OF tt-KINKS 



11 



It follows that M^) = Oij-") for some tr > if = C'(7"'^), i.e., if x = a\ In 7] at 
leading order. Using this, it is a straightforward procedure to show that the rest term 
j^R2{x;j) in (3.12) is of 0(72-2'^) for a; = a| ln7| + 0(1) (and a > 0). Hence, the 
approximation of (phix) by expansion (3.12) breaks down as x becomes of the order 
I ln7|. On the other hand, it also follows that (j^lpp^ix) = (j)fiux{x)+j4'i{x) is a uniform 
0(7)-accurate approximation of (phix) on an interval (—00, L] for L = i| ln7| + 0(1). 
Since </)flux(-L) + l4'iiL) = 0{y/j) for such L, we can compute = ^\ \n"f\ + 0(1), 
as L-j^ is the value of x at which 

= ^Mx) = ^4>lppAx) + 0(7) ^ ^0fiux(a;) + l^cpi\u{x) + 0(7). 
ax ax ax ax 

We introduce Y hy = so that it follows by (3.2) and (3.17) that F = + 
0(V7), i-e. 

i47) = i|ln7|+ln^ + 0(V7)- 

A straightforward calculation shows that (in the new coordinates) 

(t)h{L^) = 2tt- 2V¥V7 + 0(7). 

As (j)h{x) and 4>fim(x) both converge exponentially fast to fixed points which are 
order 7 apart for x —00, it follows immediately that (?ii(a;; 7) = 0(1) and R2{x; 7) = 
0(1), uniform for .t < 0. □ 

Now we are ready to consider the stability of the various types of 7r-fluxons indi- 
vidually. 

3.1. Stability of the type 1 solution. Lemma 3.5. For all < -f < jcr, 

all eigenvalues of £^{x;^) are strictly negative. For 7 = 7cr, the operator £^{x]jcr) 
has as its largest eigenvalue. For 7 = 0, the largest eigenvalue is 1). 
Furthermore, for all < j < jci, the type 1 semikinks (j)l.{x;j) are Lyapunov stable 
in the following sense. For all e > 0, there is some 5 > Q such that any solution 
(l){x, t) of the semi-fluxon equation (3.1), which is convergent to at x ^ —00 and to 
TT at X ^ +00 and which satisfies initially \\4>i'jO) ~ 4>'!r{'t'~f)\\Hi + \\(j>t{',0)\\L2 < ^ w^^' 
satisfy U{;t) - <PU-,j)\\l, + ||<^t(-,i)lk. < e for all t € M. 

Proof. Prom Lemma 3.3 it follows that has an eigenvalue A = at 7 = 7cr- The 
eigenfunction is Dx'^i(2;; 7cr) and this function is always positive, since (/)^(x;7cr) is 
monotonically increasing. From Sturm's Theorem, it follows that A = is the largest 
eigenvalue of £^ at 7 = 7cr- Next we consider 7 = 0. We can explicitly determine all 
eigenvalues of jC^{x;Q). From the explicit expression for 0^ it follows that £^(x;0) is 
a continuous even operator. For fixed A, the operator £^(a;; 0) — A has two linearly 
independent solutions. Since the fixed point is a saddle point and the decay rate to 
this fixed point is like e~^, there is one solution that is exponentially decaying at 
+00 and there is one solution that is exponentially decaying at —00, if A > — 1. If we 
denote the exponentially decaying function at —00 by w_(x; A), then the exponentially 
decaying function at +00 up to a constant is given by v+{x;A) = V-{—x;A) (since 

is symmetric in x). Obviously, w+(0;A) = V-{0]A), hence A is an eigenvalue if 
D^v+{0;A) = D^v^{0;A), (i.e., when D^v_{0;A) = 0) or if w_(0;A) = 0. 

Using [22], we can derive explicit expression for the solutions V-{x; A) (see also [10]). 
Using xi = ln(-\/2 + 1), we have 



t;_(x;0) = sech(x - a;i), u_ (x; A) = e'^'^^ [tanh(a; - xi) - yu], /U = VA + 1. 



12 



BERKS ET AL. 



Fig. 3.2. (a) The eigenvalue of linear operator associated to the type 1 semifluxon as a function 
of the bias current^. The dashed line is the boundary of the continuous spectrum, (b) A simulation 
of the evolution of a w-kink in the presence of a bias current above the critical value (7 > -fci)- The 
plot presents the magnetic field 4>x ■ The numerics show that the instability leads to the release of 
wave trains of traveling wave fluxons. In this evolution a damping, which is proportional to cpt, has 
been applied to the system. 



A straightforward calculation shows that v- (0; A) ^ 0. The condition D^V- (0; A) = 
gives that 

11 1 1 

Lt^ - -V2u- - =0, hence VA + 1 = - V2(V5 - 1) ^ A = — (V5+1). 
2 2 4 4 

Now assume that the operator C^{x;j) has a positive eigenvalue A^(7) for some 
< 7 < 7cr- Since A depends continuously on 7, there has to be some < 7 < 7cr 
such that A^(7) = 0. However, from Lemma 3.3 it follows that this is not possible. 

Nonlinear or Lyapunov stability can be derived by looking at the "temporal 
Hamiltonian" 



\p^ + \{cj>.f-cos{(j> + e)-^{ci, + e) 



dx 



This fimctional is a Lyapunov function for the system (3.1), i.e., any solution t) G 
if2(M) of (3.1) satisfies j;{H{(t>,(l>t) = 0, hence W((/)(-, t), .^((^ t)) = H(0(-, 0), 0f (•, 0)) 
for any < G M. Furthermore, the linearization D'^H at {(j>,p) = {(j)]^,^) (the point 
related to the 7r-fluxon) is given by 



-c\x-n) 
I 



which is a strictly positive definite self-adjoint operator on L2(R) x L2(R) with domain 
if2(K) X i2(R). So there is some c > such that for any {(jj^p) G i72 x L2, we 
have H{^,p) - Hi^l,0) > c(||</. - </.i + \\p\\lj, see e.g. [15, 36]. Finally, it is 
straightforward to prove that there is some C > such that H{(j),p) — iJ((/)^,0) < 
^^(110 - <f>l\\k + IbllD for any (0,p) e H2 x L^. □ 

Using standard procedures in MATLAB, the eigenvalues of the type 1 7r-fluxon 
have been calculated numerically as a function of the applied bias current 7 and 
are presented in Figure 3.2(a). Further details of the computational procedure are 
presented in section 6. Figure 3.2(a) shows that the type 1 semifluxon has only one 
eigenvalue. This eigenvalue tends to zero when the bias current 7 approaches the 
critical value 7cr as has been derived analytically. It was first proposed in [16, 20] 
that a constant driving force can excite the largest eigenvalue of a semifluxon toward 
zero. 

When we apply a bias current above the critical value 7cr, numerics show that the 
stationary 7r-kink bifurcates into a semifluxon that reverses its polarity and releases a 
fluxon. This process keeps repeating itself: the semifluxon changes its direction back 
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and forth with releasing a fluxon or antifluxon in every change. A simulation of the 
release of fluxons from a scmifluxon is presented in Figme 3.2(b). In experiments, the 
polarity of a scmifluxon can also be reversed by applying a magnetic field [14]. 

When 7 = 7(,j., the type 1 and type 3 semifluxons coincide. Prom the numerical 
analysis of the eigenvalues of the type 3 semifluxon (sec section 3.3 for details), it 
follows that there is an eigenvalue at the edge of the continuous spectrum for 7 = 7(.r. 
We conjecture that this eigenvalue bifurcates into the edge of the continuous spectrum 
at this point as 7 increases to 7cr (see Figure 3.4). 

3.2. Instability of type 2 solutions. Lemma 3.6. For all < j < 7*, the 

largest eigenvalue of C^{x;^) is strictly positive. In the limit 7^0, the largest 
eigenvalue ofC'^{x;j) converges to j{\/5 — 1). 

Proof. Using the approximation for the homoclinic orbit (f>h{x;"f) in Lemma 3.4, 
we see that, for 7 small, an approximation for the 7r-fluxon of type 2 is given by (as 
before, xi = ln(l + \/2)) 

f ^ifluxl-T + xi) +0(7], ^ x<0 

(f)l{x;'y) = < n + (pau^{x) +j(j)i{x) +j^R2{x;j), < x < L^{"/) + xi 

[ 7r + ^flux(-^) +7</'i(-^) +7^-R2(-S;7), x>Lj,{j)+Xi 

(3.18)^ 

with X = X — xi and x = x — 2L^{j) — xi. 

There is no limit for 7 ^ 0, since the semifluxon breaks in two parts, one of them 
being the 37r-fluxon 0§7r(a;; 0). In a similar way as we found the largest eigenvalue for 
the linearization operator C^(x; 0) about the 7r-fluxon (f)]^{x; 0), we can find the largest 
eigenvalue for the linearization operator £^(a;;0) about the 37r-fluxon 03^(a;;O). The 
largest eigenvalue is A^(0) = j{^/5 — 1) and the eigenfunction is 



ij^{x;0) 



e^»("^+^i)(Aio-tanh(a; + a;i)), x < 0, 

e''o(-^+^i)(/io -tanh(-a; + a;i)), a; > 0, 



where iio = \/A2(0) + 1 = ^^2(1 + ^5). (It can be shown that there is another 
smaller eigenvalue A = — ^ and similar eigenfunction if = = tanh(a;i), see 

Remark 3.8.) 

In a similar way, using the approximation (3.18) for (j)'^{x;j), the eigenfunction 
of an eigenvalue of 0^ for 7 small is approximated by 

e''(^+^i)(At-tanh(a; + a;i)) +0(^), a; < 

i;'\x;^) = \ fc2e-^^(/i-tanh(-^)) + /c3e''^(M-tanhJ) + 0(^), < x < L^{-f) + xi 
fc4e^(-®)(Ai + tanhx) +0(77), x > L^i-f) + xi. 



where ki and fi have to be determined. The eigenvalue A follows from ^ = VA^ + 1. 
Note that the secular term which is growing at infinity with the multiplication factor ^3 
is included in this approximation. When 7 = and fcs = 0, the first two lines in the 
definition of tp"^ are the eigenfunction of the linearized problem about the heteroclinic 
connection between and 3tt, as presented above. When 7 is nonzero, ^3 can be of 
order 0{j'^) for a > j as the secular term is of order 0(7"^/^) at a: = L,r(7) + xi. 

The constants /e2, and ^4 and the parameter fj, have to be chosen such that 
for 7 > (but small) the function ^p^{x,j) is continuously differentiable at a; = and 
x = L^{'y) + xi. Prom the continuity conditions at a; = 0, we obtain: 



4/x(/x-l)(/x+l) 
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Fig. 3.3. (a) The eigenvalues of the linear operator associated to the type 2 sernifluxon as a 
function of the bias current 7. The result that the largest eigenvalue is always positive shows the 
instability of type 2 sernifluxon. When 7 ^ 0, A ^ i(i/5 — 1) which is the largest eigenvalue of a 
3iT-kink, At ■y = 0, one eigenvalues comes out of the edge of the continuous spectrum (dashed line), 
(b ) The evolution of a Sn-kink (3. 11) for 7 = 0. The plot is presented in terms of the magnetic field 
(px ■ The separation of a fluxon from the semifluxon can be clearly seen. 

_ (3 + 2^/2)^(2M^-M^/2-l)(2M-^/2) 
4m(m^ - 1) + 

From one of the continuity conditions at a; = (7) + Xi, we obtain ki = 
^4(^2, fcs, /Lt). Now we are left with one more matching condition. Values of /z for 
which this condition is satisfied correspond to the eigenvalues of the operator C^{x; 7) 
for 7 small. More explicitly, the spectral parameter has to satisfy the equation 

(3.19) jr(^) = 16''/e3(/x - l)2(77r)-''((3/x + 4)7r7 + IS/i) + 0{^-^'+^) = 0. 

Note that this expression is not defined at 7 = 0. This corresponds to the singularities 
in the expression for (ji)^ as 7 ^ due to the fact that ^^(7) 00 for 7 — » 0. 
Evaluating ^(/u) (77r)'^ at 7 = 0, we see that there are four positive roots for fi, 
leading to four squared eigenvalues, namely A(0) = |(a/5 — 1), — i, and the double 
eigenvalue A(0) = 0. The first two come from the zeros of kz and are related to the 
eigenvalues of the 37r-fluxon. The double zero eigenvalues are the eigenvalues of the 
fluxon. One can also notice that there is no term with a multiplication factor k2 to 
this leading order. This term appears at most of order 0(7^+^). Finally, as with the 
type 1 semi-fluxon, the root fx = corresponds to the edge of the continuous spectrum 
and the "eigenfunction" is not in H2{M.). 

The proof that the largest eigenvalue is near j (a/5— 1) for 7 small will be complete 
if we can show that J^n{V^/A{l + ^ 0, i.e. the non-degeneracy condition that 
says that the eigenvalue can be continued continuously for 7 small. 

Simple algebraic calculations give that 

(3.20) + ^)) = ci7-^''+^) + 0(y-^(^+^)) 

with ci a positive constant. Hence, jr^(^(l + a/S)) > 0. 

This completes the proof that the largest eigenvalue is near ^(Vs — 1) for small 
but positive 7. Since the largest eigenvalue depends continuously on 7, it can only 
disappear at a bifurcation point. There are no bifurcation points and it is not possible 
that the eigenvalue becomes (sec Lemma 3.3), hence the largest eigenvalue will be 
positive as long as fluxon (/)J(a;;7) exists, i.e., for < 7 < 7*. □ 

REMARK 3.7. We cannot use a comparison theorem, because (fy^ < (f)^ for x <0 
and (j)^ > for x > 0. 

To consider the relation between the eigenvalues of £?{x]^) and the stability 
problem of 4>'^{x\ 7), we denote the largest eigenvalue of £^(a;;7) by A^(7). The asso- 
ciated eigenvalues for the linearizations are solution of the equation — A^(7) = 0, 
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hence A = ±yjK^{"f). Since A^(7) > 0, this imphes that one of the two eigenvalues 
has positive real part, hence the 7r-fiuxons of type 2 are unstable. The numerically 
obtained eigenvalues of semifluxons of this type as a function of 7 are shown in Fig- 
ure 3.3(a). In the proof of Lemma 3.6 we have found three different eigenvalues for 
7 small and the possibility of a fourth eigenvalue coming out of the continuous spec- 
trum at 7 = 0. In Figure 3.3(a), we see the continuation of those eigenvalues. In 
Figure 3.3(b), we present the evolution of a 37r-kink (3.11) which is the limit of a type 
2 scmifluxon when 7^0. The sciparation of a fluxon from the semifluxon is clearly 
seen and indicates the instability of the state (which confirms the analysis in the proof 
of Lemma 3.6). 

REMARK 3.8. A type 2 semifluxon can be seen as a concatenation of a 37r- and a 
— 27r-kink in the limit 7^0. In that limit the other eigenvalues of C'^{x; 7) converge 
to 0, — ^, and —1. The eigenvalues and —1 are contributions of the — 27r-kink. The 
eigenvalue — ^ corresponds to the first excited state of the 37r-kink with eigenfunction 



where /U = \/A + 1 = -4 



e'^(''^+^i)(Ai-tanh(a; + a;i)), a; < 0, 

e''(-^+^i)(tanh(-a;-|-xi)-/x), a; > 0, 



3.3. Instability of type 3 solutions. Lemma 3.9. For all < -f < 7cr, the 
largest eigenvalue of C^{x;^) is strictly positive. For 'y = 7cr, the operator C^{x;jcr) 
has as its largest eigenvalue. 

Proof. The solution (/)^(x,7cr) = 0^(a;,7cr), hence from Lemma 3.5 it follows that 
the largest eigenvalue is A = 0. 

For 7 near zero, we will use the approximation for the homoclinic orbit (j)h{x;^) 
in Lemma 3.4 to get an approximation for the type 3 fluxon 



Z^'appr (S ) = 'jf'flux (^ ) + T^^i (^ ) + 7^ -R2 (S ; 7) , a; < - L„ (7) + a; 1 , 

f^lpprix) = (t>Ru4-x) + iM-x) + ■y'^R2{-x; 7), -Kij) + a;i < a; < 0, 
^lppr{-x - a;i) = TT + ^fiux{-x - xi) + 0(7), a; > 0, 



where a; = a; — xi and x = x — xi + 2X7^(7). 

In the limit 7^0, the type 3 semi-fluxon break into a type 1 semi-fluxon and 
a fluxon. Both are stable and the largest eigenvalue of the fluxon is zero, while the 
largest eigenvalue of the type 1 semi-fluxon is negative. Hence to approximate the 
largest eigenvalue of the type 3 semi-fluxon for 7 small, we set 

A(7) =7Ai(7). 

To construct the first part of the approximation of the eigenfunction, we consider 
X < —L^(pi) +Xi, i.e., X < Ltj(pi). In this part of the arguments, we will drop the hat 
in X. On (—00, Z/jr), we expand V'ipprox = i^o + li^i-, this yields the following equations 
for Vo,i(a;), 

(3.21) £Vo = 0, £Vi = [Ai(0) - .^i(a;)sin,^flux(a;)]Vo. 

As V'ipprox to be an eigenfunction, we have V'approx(^) as x ^ — 00. Further- 
more, we remove the scaling invariance by assuming that V'approx(O) = 1- This implies 
that iPq{x) is given by 
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(see 3.16). To solve the ^Z)! -equation, we note that j^(j)i{x) is a solution of 

C-i) = -(j>i sin^flux^?!'flux = -2?!>i sm.(j)^^^tj)o, 
(see (3.15) and (3.2)) so that we find as general solution, 

^.(.) = - iA.(ln(cosh.) + £ 

( \ + sinhx) + 
V cosh X I 2dx 



+ 



B + ^Ai tanha; 



Using lima,_^_oo '0i(a;) = and V^ilO) = we find that A = \, B = |Ai(0). As in 
the case of 0i (x) , we are especially interested in the unbounded parts of ■i/'i (a;) and 



(3.23) 



i\}y\u(x) = iAi(l + tanhx) sinha; — ^ arctane"^ cosha;, 
^V'i|«(a;) = ^Ai(l + tanhx) coshx — I arctane^ sinha;. 



We note that the error term \tp{x) — i^lppiix)] = 7^|S'2(a;; 7)| is at most 0(7) on 
(— (the analysis is similar to that for 7^|i?2(a;; 7)]). 

Next consider the second part of the approximation, i.e., x between —L^{j) + 
Xi and 0. Here we define the translated coordinate x = x — Xi, which is on the 
interval (— Ltt, — a;i) and again we drop the tildes. Since we have to match V'appr(a^) to 
the approximation ip'^ppj.[x) of ^{x), along (/)^ppj.(x) and thus defined on the interval 

(— Ltt,— xi), we need to compute -tplppTi^-^) ^^'^ ^^appr(-^7r) which to the leading 
order are calculated from (3.23), i.e. 

(3.24) i'lppAL.)='-^V^+0{j), i^l..iL.)= ^^'^^~" Vl + 0{^). 

Thus, both Vippr(i^) and £V'appr(i7r) are 0(^7)- 

Now, we choose a special form for i/'^pp^(x), the continuation of 'ip{x), i.e. the part 
linearized along 0^ppj.(a;). It is our aim to determine the value of Ai, for wliic;li there 
exists a positive integrable solution tjj of C'\x; 7)?/) = ^Ai{Q)'if'j. By general Sturm- 
Liouville theory [33] we know that this value of Ai must be the largest eigenvalue. Our 
strategy is to try to continue ip{x) beyond (—00, Lj^) by a function that remains at 
most 0{y^), i.e. we do not follow the approach of the existence analysis and thus do 
not reflect and translate '^i'appr(2^) to construct V'apprC^;) (since this solution becomes 
in general 0(1) for x — 0{1)). Instead, we scale ip'^pp^ix) as ■yip{x). The linearization 

tl){x) along (/>appr(2;) on the interval {—Lt^jXi) must solve £tp — 0{-y), thus, at leading 
order 

(3.25) i(i(x) = + B (—^ + sinhx) . 

cosh x \ cosh x I 

The approximation V'appr(a^) = li'i'A must be matched to ■i/'ippr(-^7r) and ^i/)^ppj.(Lt) 
at a; = — L^, i.e. 
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Note that A does not appear in these equations; as a consequence, tplpp^ix) and 
tp'ipprix) can only be matched for a special value of Ai, Ai(0) = jtt, with B = 
— Ai(0) < 0. Thus for this special vahic of Ai and for A > 0, wc have found a positive 
C^-continuation of the solution 'iJj{x) of the eigenvalue problem for C^{x;^) - recall 
that a; < in the domain of tpix). At the point of discontinuity {—xi for ip{x), or at 
a; = in the original coordinates of (3.1)), we have 

(3 26) ^'pp'^^^i) = iH-xi) = l[kV2A - |^/2(ln(y2 - 1) - ^2)] + Oij^), 
£V'fppr(-^i) = 7£V^(-xi) = 7[^i - f (ln(y2 - 1) + 3^/2)] + 0{j^). 

Hence, we have constructed for a special choice of A, A = A* = f 7 + C(7y^) > 0, an 
approximation of a family of positive solutions of the eigenvalue problem for C^{x; 7) 
on a; < in the coordinates of (3.1) - that attain the values given by (3.26) at 
a; = 0, and that decay to as a; ^ —00. The question is now whether we can 'glue' 
an element of this family in a -fashion to a solution of the eigenvalue problem for 
C^{x;'y) on a; > - with A = A* that decays (exponentially) as x 00. If that is 
possible, we have constructed a positive integrable solution to the eigenvalue problem 
for £^(a;;7), which implies that A* > is the critical eigenvalue and that ?!>^(a;) is 
unstable. 

An approximation of il^{x) on x > 0, ^'Ipprix), is obtained by linearizing along 
(jf)^ppr(a;) and by translating x so that x G {xi,oo). Since V'appr(^) has to match to 
expressions of 0(7) (3.26) at Xi, we also scale V'appr(^)) ''Ptppiix) = 7'4'{x)- We find 
that Cip = 0(7) so that ■i/j{x) again has to be (at leading order) a linear combination 
of i^bix) and ipuix) (3.16). However, ip must be bounded as a; ^ 00, which yields that 
V'(a;) = A/ cosh a; + 0{j) for some A e ffi. At the point of discontinuity we thus have 

,3 27x i^lpprixi) = liKxi) = ^V2Aj + 0{j^), 

l^lppAxi) = l^M^i) = -ii7 + 0(7^). 

A positive C^-solution of the eigenvalue problem for £^(a;;7) exists (for A = A*) if 
there exist A,A>Q such that (sec (3.26) and (3.27)) 

iV2i - f^/2(ln(^/2-l)-^/2) = iy^i 
\A - f(ln(V2-l) + 3^/2) = -\A 

Since the solution of this system is given hy A = jTt[\/2 + ln(\/2 — 1)] > and 
A = i7ry2 > 0, we conclude that the eigenvalue problem for the 7r-fluxon 4>^{x;^) 
has a positive largest eigenvalue 

(3.29) A, = ^7 + 0(7/7). 

Hence the eigenvalue for 7 small is positive. From Lemma 3.3 it follows that there 
are no zero eigenvalues between and 7cr, hence the largest eigenvalue of £"^(7) is 
positive for all values of 7. □ 

REMARK 3.10. For any A = 0{^), or equivalently any Ai = 0(1), there exists 
a (normalized) solution to the eigenvalue problem for £^{x; 7) on a; < that decays 
as X ~oo, and that is approximated by 4'lppiix) and '0^ppj.(x) (matched in a C^- 
fashion at ±L„). If Ai is not 0{^) close to jw, however, '0appr(^) cannot be scaled 
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Flc. 3.4. The eigenvalues of the the linear operator assoeiated to the type 3 semifluxon as a 
function of the bias current 7. The result that the largest eigenvalue is always positive shows the 
instability of type 3 semifluxon. When 7 <g; 1, according to (3.29) the largest eigenvalue is approx- 
imated by A = ^7 shown in dash-dotted line. The dashed line is the boundary of the continuous 
spectrum. 

as ji^ix) and the solution is not 0(7) at the point of discontinuity in general it is 
0(1). Moreover, for any Ai = 0(1), there also exists on a; > a l-parameter family of 
(non-normalized) eigenfunctions for the eigenvalue problem for C^{x; 7) that decay as 
a; ^ 00. In this family there is one unique solution that connects continuously to the 
(normalized) solution at a; < 0. In fact, one could define the jump in the derivative at 
X = 0, J'(A;7), as an Evans function expression (note that J{\]^) can be computed 
explicitly at 7 = 0, see [10]). By definition, is an eigcnfunction of £^(a;;7) if and 
only if i7(A;7) = 0. In the above analysis we have shown that J'(A*;7) = for 

REMARK 3.11. The classical, driven, sine-Gordon equation, i.e. 9 = and 7^0 
in (3.1), has a standing pulse solution, that can be seen, especially for < 7 -C 1, 
as a fluxon/anti-fluxon pair. This solution is approximated for -^(j) > (the fluxon) 
by (/>appr(2^) and for ^0 < (the anti-fluxon) by (/)appr(— a;). It is (of course) un- 
stable, the (approximation of the) critical unstable eigenvalue can be obtained from 
(3.24). The corresponding eigcnfunction is approximated by V'appr(2;) o'^ (— oo,^^), 
and we conclude from (3.24) that |^Vippr(-^7r) = for A^ = 7A1 = 7f + 0(7^7) 
(while 7/;;|^pp^(L^) > 0). Hence, for this value of Ai, we can match V-'appr(-^) to 
V'appr(^) = V'appr(~^) ^ C^-fashion, it gives a uniform C'(7)-approximation of 
the critical, positive (even, 'two-hump') eigenfunction of the fluxon/anti-fluxon pair 
at the eigenvalue A+ = ^\/27r^/7 + 0(7) > 0. 

To consider the relation between the eigenvalues of £^(a;;7) and the stability 
problem of 4>%{x;'y), we denote the largest eigenvalue of C^{x\^) by A^(7). The 
associated eigenvalues for the linearizations are solution of the equation A^— A'''(7) = 0, 
hence A = ±\Jts?(j). Since A^(7) > 0, this implies that one of the two eigenvalues 
has positive real part, hence the fluxons of type 3 are unstable. In Fig. 3.4, we present 
numerical calculations of the eigenvalues of type 3 semifluxon as a function of the bias 
current 7. 

REMARK 3.12. A type 3 semifluxon can be seen as a concatenation of a 27r- and 

a — TT-kink in the limit 7^0. In that limit the other eigenvalue of C^{x\ 7) converges 
to -i(V5+ 1) (Figure 3 .4) which is a contribution of the — 7r-kink. 

4. Lattice 7r-kinks and their spectra in the continuum limit. In this 
section, we consider (2.9) for a small lattice spacing a, i.e., the driven O-tt sine-Gordon 
equation with a small perturbation due to lattice spacing effects. For a = 0, the 

semifluxons of all types are constructed as heteroclinic connections with transversal 
intersections at .x = in the two-dimensional phase space of the static equation (2.10). 
Therefore, all three types of semifiuxons will still exist in the perturbed system with 
< a <C 1, see [11]. The three types of semifluxons are denoted as (j)\{x;a;^), for 
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Fig. 4.1. The phase portrait of the stationary system (2.10) for 7 = and some values of the 
lattice spacing a. The dashed lines are the unperturbed phase portrait for o = and the other lines 
correspond to a = 0.5. 



i = 1, 2 and 3. In Fig. 4.1, we present the phase portraits of the sine-Gordon equation 
both with and without the effect of a perturbation due to lattice spacing. 

The lattice spacing a does not affect the stationary points of the phase portraits, 
as can be easily checked. The existence parameters 7* and 7cr will be influenced by 
the lattice spacing a. For a small, they are 

^ V ' 

!»0.0109 

, , , , 2 Vtt^ - 4 - tt + 2 arcsin(^) „ , 

^ ^ 

^0.0223 

see [31] for details. For 7 > 7cr(a) no static semifluxon exists. 

As we have seen in the last section, for o = 0, the type 3 semifluxon is marginally 
unstable at 7 = 7cr and 7 near zero. So there is a possibility that lattice spacing effects 
stabilize the type 3 semifluxon near those values of 7. However, it turns out that this 
is not the case and the stability of the semifluxons is similar to the case a = 0. 

Theorem 4.1. For a small, the linearizations about the w-kinks have the following 
properties: 

(i) The eigenvalues of the linearization about the monotonia type 1 n-kink (/)^ (x; a; 7) 
are strictly negative for < 7 < 7cr(a). At = 7cr(a); the largest eigenvalue is zero. 
These n-kinks are linearly stable. 

(ii) The largest eigenvalue of the linearization about the monotonia type 2 n-kink 
(f>'^{x; a; 7) is striatly positive for < 7 < 7*(a). These n-kinks are linearly unstable. 

(iii) The largest eigenvalue of the linearization about the monotonia type 3 n-kink 
(j)^{x; a; 7) is strictly positive for < 7 < 7cr(«)- These n-kinks are linearly unstable. 
In the limit for 7^0 and 7 7cr(a), the largest eigenvalue converges to zero. 

The proof of this theorem will proceed along similar lines as the proof in the 
previous section. First we consider the eigenvalue problem of a solution (j))^(x]a\^), 
which can be written as 

a; 7) w = w, 

where D{x; a; 7) is now defined as the linearization associated with (2.10), i.e., 

C'{x;a;^) = Da:x - cos {(pl{x-,a;j) + e{x)) 

- 2cos(^i:>^^ - 2((/);(a;;0;7)):rSin(^i:>^ 

-1 + 27 sin - {{^l{x; 0; 7)),)^ cos + 0{a^) 



20 



BERKS ET AL. 



where (j) = 4'\r{x] 0; 7) + 9{x). 

Lemma 3.3 can be extended to a 7^ and give a necessary and sufficient condition 
for C'{x; a; 7) to have an eigenvalue A = 0. 

Lemma 4.2. The eigenvalue problem 



has an eigenvalue A = if and only if one of the following two conditions holds 

(i) Dxx<l>lr{x'i^'^) is continuous at X = 0, i.e., <j)l^{0;a;^) = kn — a^j^ ^{(^^)> 
for some k G Z; 

(ii) Dx(j)l.{0;a;^) = and there are some x±, with sgn(a;±) = ±1, such that 



Proof. As the proof of Lemma 3.3 is based on the fact that the derivative of the 
semifluxon is a solution of the linearized system for x 0, we can follow the same 
arguments to prove this lemma. Again this leads to two conditions that either is 
continuous at a; = or the second condition as stated above. 

In order to determine when is continuous, we use the static equation (2.10) 
and expand near a = 0: 



again with (f) = (/>^(x;0;7) + 6{x). The continuity of D^x^lr at a: = leads to the 
expression for 0^(a;;a;7) as given above. □ 

At 7 = 7cr(a), the stable manifold of the tt + arcsin7 and the homoclinic con- 
nection at arcsin7 are tangent, implying that Dxx4'\{^A0''t1) is continuous at 3; = 0. 
Thus the first condition of the lemma is satisfied at 7 = 7cr(«) for i = 1,3. For the 
same reasons as before, the second condition is never satisfied. 

Since A = is an eigenvalue of the linearized operator £*(x;a;7) if and only if 
7 = 7cr(a), the sign of the eigenvalues of £*(a;; a; 7) will not change. Thus the behavior 
of the eigenvalues near 7 = will again determine the stability of the semifluxons. 

For 7 = and 6 = 0, the sine-Gordon equation with a perturbation due to the 
lattice spacing has a heteroclinic orbit connecting and 2it. As before, the heteroclinic 
orbit will play an important role in determining the stability of the semifiuxons for 
small values of 7. For small values of the lattice spacing a, we can approximate this 
heteroclinic orbit up to order by using the 27r-fluxon (?iiflux and its linearization. 

Lemma 4.3. Let (^§^^(0;) denote the heteroclinic orbit of the sine-Gordon equation 
with a perturbation due to the lattice spacing (i.e., (2.9) with 6 = and 7 = 0^. For 
the lattice spacing a small, we have for the symmetric (i.e., 0g^^(O) = n) heteroclinic 
connection (f)^^^{x) 



C{x;j)v = Av, xeR, 



Dx<l>i{x±;an)^0. 




(4.3) 




(/)fln^{x) + a^<l)a{x) + 0{a'^), 



where 



(4.4) 



1 — 3sinha; + a;cosha; 



12 cosh^ X 



This approximation is valid, uniform in x €M. 
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Proof. The spatially localized correction to the kink shape ^flux(;c) due to the 
perturbation term representing lattice spacing is sought in the form of perturbation 
series: 

'/'fluxCa;) = 0flux(x) + aVa(x) + 0{a'^). 
It is a direct consequence that (j}a{x) satisfies 

{^.bp-{x;Q)(l}a{x) = f{x) = -j^ [2cos(/)flux(a;)92:x<^flux(a^) 

- sin^flux(a;)(9a;(/)flux(a;))^ - cos(/)flux(a;) sin^flux(a;)] , 

where £^{x;0) is the linearized operator associated to the fluxon, i.e., jC^{x;0) = 
Dxx - cos0flux(a::)- 

Using the variation of constants method, we obtain the general solution of (4.5), 

i.e. 



(4.6) 
where 

A{x)- 
B{x)- 



(l>a{x) = A{x) secha; + B{x) (a;secha; + sinha;), 



Ao + 



1 

24 
1 

24 



2 In 
2 + 



1 — cosh x — sinh x 



\ cosh X 
1 

cosh^ X 



1 — sinha; 
3 

cosh* X 



+ 



6 sinh X 4 sinh x 



cosh X cosh^ X 



f 

Jo 



mi 

cosh^ 



The integration constant Bq is determined by the condition that (paix) is bounded, 
leading to -Bq = j2- The integration constant Aq is determined by the requirement 
that 0flux(O) = hence ^a(O) = 0, giving that Aq =0.D 

For 7 = 0, the static model (2.10) for a O-tt Josephson junction with lattice spacing 
effects has both a tt- and a 37r-kink solution. The 27r-heteroclinic orbit found above, 
can be used to derive approximations for those kinks. 

Lemma 4.4. For a small and j = 0, we have an explicit expression for the n- 
and Sn-fluxon up to order 0{a^), respectively: 



(4.7) 



-ul{x - ln(l + V2)), for 

774(-.T-ln(l + V2)), for 

-ul^{x + \n{l + V2)), for 

ul^{-x + \n{l + V2)), for 



x<Q 
x>Q 
a; < 
a; > 



0i(x;O) + a2 
4>l^{x;a;Q) = <^i„(a;; 0) + 

where </>i(x; 0) and (pl^^ix; 0) are the tt- resp. the Stt- fluxons as defined in (3.11) and 
= 12^ (¥ - ^ V2) + 3 tanha; - x) 

(-^ + 5 ln(3 - V2) + 3 tanha: - a;) . 



12 cosh X 



4.1. Stability of type 1 semifluxon. We will show that the type 1 wave 
(/)^(a;;a;7) is linearly stable for small a and < 7 < 7cr by analyzing the largest 
eigenvalue of jC^{x; a; 7) for < 7 < 7cr(a)- 

Lemma 4.5. For the lattice spacing parameter a .sufficiently small and < 7 < 
7cr(a), the largest eigenvalue ofC^{x;a;j) is strictly negative. For ^ — 7cr(a), the op- 
erator {x; a;7cr(a)) has as its largest eigenvalue. For^ = 0, the largest eigenvalue 
decreases as a increases and is proportional to — ^(v^ + 1) — 0.0652a^ + 0{a'^). 
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Proof. First we look at the stability of the 7r-kink at 7 = 0. Writing v{x) = 
v^{x) + a^i/(.-r) + 0{a'^) and A = Aq + a^Ai + 0{a!^) and expanding the eigenvalue 
problem for the stability of the 7r-kink (j)\{x;a;Q) in a Taylor series, result in the 
following equations 

(£i(x;0;0)-Ao)^;O(x) = 0, 
^^■''> {c^{x;Q;Q)-Ko)v\x) = {A, - ui{x) sm{cPi{x;0) + 9)) v^{x) - g{x), 

where /x = VAq + 1, Aq = -1(^5+ 1), 

gA'(x-in(i+\/2))[^j^^j^(^_l^(^_^y2))-/z], for a;<0, 

gP(-x-in(i+y2)) [tanh(-a; - ln(l + %/2)) - ^] , for a; > 0, 

2^0^ Ao + w° + 2wO^ cos ^(x) - 2 cos^ (^(.t)?;" - 29^,^ ((^^ (x; 0)) sin 

- 2d^(l)l{x;0)sm^{x)v° - (^^^^(a;; 0))^ cos0(x)t;O - ?;0A2 - 2uOAo cos^(x)] , 

with again (f){x) = (p^{x; 0) + ^(a;) (see Lemma 3.5). 

The parameter value of Ai is calculated by solving (4.8) for a bounded and de- 
caying solution v^{x). The general solution can be derived by using the variation of 
constant method because we have the homogeneous solutions of the equation. One 
can also use the Fredholm theorem (sec, e.g., [30]), i.e. the sufficient and necessary 
condition for (4.8) to have a solution G -ff2(]R) is that the inhomogeneity is per- 
pendicular to the null space of the self-adjoint operator of jC^{x; 0; 0). If ( , ) denotes 
an inner product in H2{M.), then this condition gives 

= ((£i (a;; 0; 0) - Ao)v\v°) = {A^v^ - sin(0i (a;; 0) + ^) - g, v"") 

which implies that 

(4.9) A, = 3584(70^^1 + 75)- 99(1 + V5)) ^ 

24576(-70\/T0 - 350V2 + 495 + 99V5) 

Now assume that the operator £^{x;^) has a positive eigenvalue A-^(7) for some 
< 7 < 7cr(a)- Since A depends continuously on 7, there has to be some < 7 < 
7cr(a) such that A^(7) = 0. However, from Lemma 4.2 it follows that this is not 
possible. □ 

4.2. Instability of type 2 semifluxon. In Lemma 3.6 we have seen that for 
a — 0, the linearization about the type 2 semifluxon has a strictly positive largest 
eigenvalue. Also the limits of this eigenvalue for 7^0 and 7^7* are still strictly 
positive. Thus a small perturbation associated with the lattice spacing can not sta- 
bilize the type 2 semifluxons. 

For completeness, we will consider the case 7 = 0. In this limit, the type 2 
semifluxon can be seen as a concatenation of a 37r-kink and a — 27r-kink. As before, 
the limit of the largest eigenvalue for 7^0 will be equal to the largest eigenvalue of 
the 37r-kink. We have seen that the largest eigenvalue of the 37r-kink at 7 = and 
o = is strictly positive and the following lemma shows that small lattice spacing 
effects increase this eigenvalue. 

Lemma 4.6. For the lattice spacing parameter a sufficiently small, the largest 
eigenvalue of the linearization £^(a:;o;0) about the Sn-kink (j)'^^{x;a;0) is strictly 
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positive. Moreover, it increases as a increases and is proportional to ^(\/5 — 1) + 
0.0652a2 + O(a4). 

Proof. Note that the lowest order analytic expressions for the tt- and the 37r-kink 
differs only in the sign of the 'kink-shift' (sec (4.7)). Because of this, we can follow the 
same steps as the proof of Lemma 4.5. Writing the largest eigenvalue of a 37r-kink as 
A = Ao + a^Ai + 0{a'^), with Aq = (v^ — l)/4 as has been calculated in Lemma 3.6, 
we compute Ai to be: 

(4.10) A, = 3584(665857(V5- 1) -470832V2(V5 + 1)) ^ ^ ^^^^ 

24576(3329285 - 2354160\/2 - 665857\/5 + 470832^/l0) 

□ 

Thus up to order 0{a*) the lattice spacing effects destabilize the 37r-kink. 

Because a 27r-fluxon in the 'ordinary' sine-Gordon equation can be pinned by 
lattice spacing effects, one might expect to have a stable 37r-kink in the O-tt sine- 
Gordon equation with larger lattice spacing effects. This is confirmed by numerical 
calculations in section 6, see Figure 6.7. If the 37r-kink is stable for 7 = 0, a stable 
type 2 semi-kink might exist for 7 > when the repelling force between the 37r- 
kink and the anti-fluxon is smaller than the energy to move a fluxon along lattices. 
However, in section 6 it will be shown numerically that the type 2 semikink is unstable 
for all values of the lattice spacing, see Figure 6.8(b). 

4.3. Instability of type 3 semifluxon. For 7 small or close to 7cr, it has been 
shown in Lemma 3.9 that the type 3 semifluxons are weakly unstable. This opens the 
possibility that the perturbation term representing the lattice spacing stabilizes the 
semifluxon. This is not the case however. 

Lemma 4.7. For small lattice spacing a and bias current < 7 < 7cr(a), the 
largest eigenvalue of the linearization C'^{x; a; 7) about the type 3 semifluxon (/>^(x; a; 7) 
is strictly positive. For 7 = 7cr(a); the operator £^(a;;a;7cr) has as its largest 
eigenvalue. For 7 near zero and c? = ^c? , the largest eigenvalue of C^{x;a;^) is 

A* = (f + Tio«') 7 + 0(7V7). 

Proof At 7 = 7cr, the solution (j)\{x;a;^„{a)) = (t)\{x;a;"f„{a)). Hence from 
Lemma 4.5 it follows that the largest eigenvalue of the linearization about 4)%{x; a; 7cr(ct)) 
vanishes. 

Prom Lemma 3.9 it follows that the largest eigenvalue of the linearization about 
(/>^(x; a; 7) is positive for a = and < 7 < 7cr- Thus a small perturbation can not 
change the positive sign of the largest eigenvalue if 7 is not near or 7cr. Now assume 
that a small perturbation would lead to a negative largest eigenvalue near 7 = or 
7 = 7cr. Then there has to be a zero eigenvalue near 7 = or 7 = 7cr, but this is not 
possible according to Lemma 4.5. Thus we can conclude that the largest eigenvalue 
is always positive. 

To complete the proof, we will derive the asymptotic expression of the eigenvalue 
near 7 = 0. Since both a and 7 are small, we relate those two parameters by writing 
= 76^. Now the approximation for the type 3 semifluxon can be written as 

{0flux(x)j- l4)i{x) +ja^(t>a{x) + 7^i?2(3;;7)> _ x< -L^i^y) + xi, 

(t)fLu^{-x) + l4>l{-x) + ^a?'(j)a{-x) + ^'^R2{-X\ 7), -^Trd) +Xi<X<Q, 
7r-|-(Aflux(-a;-a;i) + e>(7), a; > 0, 

where x = x — xx+ 2L^(7) and x = x — xi. It can be shown that the shift ^^(7) does 
not depend on a? in lowest order, i.e., I/7r(7) = 5I ln7| + In + 0(^7). 
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To find largest eigenvalue, we set again A^(7) = 7Ai(0) and follow the steps in 
the proof of Lemma 3.9 with some additional terms added to some expressions. 

First, we consider the part of the approximation with x < — L^(7) + xi or x < 
L„{j). As before, we drop the hat in x in this part of the arguments. On (—00, Lj^), 
the general solution of the eigenvalue problem of the order 0{'y) after expanding 

V'ipprox = ^0 + iV'l is 



ipi{x) 



— — -Ai I Incosha; + / -^5— ( 

4 2 'V Jo cosh^^ 



1 



coshx 



+ 



^Ai(O) + ^Ai tanhx 



V cosh a; / 2 \ax ax J 

[16 In 2 + e^" (32 In 2 - 295 + 60x) + 30x + 137 + 7e^ 



360(e2^ + 1)3 

- 16 ln(e2^ + l)(e2^ + 1)^ + e'*^(151 + 30x + 161n2)] . 

We note that the error term |'^(a;) — tplppii^)] = 7^|'S'2(a;; 7)] is still at most C(7) on 
(-00, L^). 

Next consider the second part of the approximation, i.e., x between — L7r(7) + Xi 
and or a? < — £^(7). Again, we drop the tilde in x. We scale tp{x) as ■yip{x). The 
linearization ipi^) along 0appr(^) the interval {—L„,—Xi) must solve £-0 = 0(7). 
Thus, at leading order 



ip(x) = — 4 1" B ( — ^ h sinhx^ . 

cosh x \ cosh X J 



The last part of the approximation of '^(x) on a: > 0, i\}\^^^(x\ is obtained by 
linearizing along </'appr(^) ^^"^ translating x so that x € (xi,oo). We also scale 
V'appr(^) = ^y^i^)- As must be bounded for x 00, it follows that tp{x) = 
A/ cosh a; + 0(7) for some AgR. 

Finally we have to connect all parts of the eigenfunction in a C^-fashion. This 
determines the values of Ai(0), A, B, and A as 

Al(0) = i7r+^a^ B = -Ai(0), i = i7r[\/2+log(x/2-l)], and A=^TrV2, 

thus Ai(0) > 0, i? < 0, A > and A > 0. And we can conclude that the eigenvalue 
problem for the 7r-fluxon (f)^{x; a^/7; 7) has a positive largest eigenvalue 

□ 

5. Semikinks in the weak-coupling limit. In this section we will consider 

the discrete O-tt sine-Gordon equation (2.3) when the lattice parameter a is large. The 
time independent version of (2.3) is well-known: when 7 = 0, it corresponds to the 
so-called Standard or Taylor-Greene-Chirikov map [9] and when 7 ^ 0, it is called the 
Josephson map [23]. Since we are interested in the case that the lattice spacing a is 
large, we introduce the coupling parameter e as e = ^ and the equation becomes 



(5.1) 



(i>n-£ - 2(^„ + ^n+i] = - sin(^„ + 6'„) + 7. 
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When there is no coupling, i.e. e = 0, it can be seen immediately that there are 
infinitely many steady state solutions: 

cos(fc„7r) arcsin7 + fc„7r, n = 0, —1, —2, . . . 

cos{knTr) arcsin7 + + l)n, n= 1, 2, 3, ... , 

where fc„ is an integer. The only monotone semi- kink is the solution with fc„ = for 
n G Z, thus it is natural to identify this semikink with the type 1 semikink. However, 
it is less clear which solution would correspond to the type 2 and type 3 semi-kinks. 
Possible candidates for the type 2 wave are solutions for which there is some N E N 
such that A;„ = for n < and n > N and A;„ = 1 for < n < A^. Similarly, 
candidates for the type 3 wave are solutions for which there is some N such that 
kn = ior n < —N and n > and kn = I ioY —N < n < 0. But there are many other 
candidates involving combinations of fc„ or fc„ = 1 as well. If one starts with 
such a wave in the uncoupled limit, i.e., with e ^ 1 or a — > oo and uses continuation 
to follow this wave in the discrete system (5.1) towards a = or £ ^ oo, then it 
turns out that most waves end in a saddle-node bifurcation. More details about the 
continuation can be found in section 6. 

In this section we will focus on the analytical study of the type 1 semi-kink for 
the coupling parameter s small (thus the lattice spacing a large). We will denote this 
wave by ^'^(n; e; 7) and for e = 0, we have 

arcsin7, n = 0, — 1,— 2, .. 



$1/^.0.^) = / arcsin7, n = 0,-l, 

*7rl".u,7; I 7r-harcsin7, n = l,2,3. 



The existence of the continuation of (5.2) for small coupling e is guaranteed by 
the following lemma. 

Lemma 5.1. The steady state solution <I>^(n;0;7), representing the semifluxon 
of type 1 in the uncoupled limit e = 0, can be continued for e small and 7 < 1 . It is 
given by 



' arcsin7 -I- ©(e ), n < —1; 



(5.2) $i(n;£;7)=< 



arcsin7 -I- e \-0(s'^), n = 0; 

TT + arcsin7 — e \-0(e^), n = 1; 



, TT + arcsin7 + O(e^), n > 2. 

For 7 close to one, we write 7 = 1 — £7. // 7 > tt, then the type 1 solution is: 

( f - Vgv ^+0(£ ), n<-l; 

(5.3) <PAn,s,l £7)- ^_^y2(^ + 0(£), n=l 

[ ^-v^v^ + 0(£), n>2 



2 

From (5.3) we obtain the critical bias current for the existence of static semifluxon as 
(5.4) 7„ = l-£7r + 0(£2). 



Proof. The existence proof for 7 < 1 follows from the implicit function theorem 
as given in [21, Theorem 2.1] or [24, Lemma 2.2]. 
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For the case 7 = 1 — £7, the impUcit function theorem as presented in the refer- 
ences above can not be appUcd immediately. However, after some manipulations, the 
impHcit function theorem can be appHed agaim First we substitute into the steady 
state equation 7 = 1 — £7 and $ = $0 + v^^, where ^oin) = f , for n < and 
$o(n) = for n > 1. This gives the following equations: 



= cos(Ve^(n))-i + ^_^[^(^_i)_2$(n) + $(n + l)] =:F„(^,£), n 7^ 0, 1 
= cos(VF^(0))-i _ 2$(o) + ^(1)] _ ^ =.^Foi^,e), n = 

= cos(Ve^(i))-i ^ ~ _ ^[^(0) _ 2^(1) + $(2)] + TT =: f 1 £), n = l 



Using that lim£_>o ■ 



s(Vi^(K))-l 



= — i($(n))^, the definitions for F can be smoothly 



extended to £ = too. The equations for £ = become 

$2(n) =27, n^^O,!; 5^(0) = 2(7 - tt); and $^(1) = 2(7 + tt). 

For |n| large, wave should be asymptotic to the center point of the temporal dynamics, 
hence $(n) = —\f^ for |n| large. So for 7 > tt, there are two monotone semi-kinks 
(recall that the full semi-kink is given by $0 + a/e^)- 



$±(n;0;7) = < 



/27, 



n < -1; 



±V ^2(7-7r) , n = 0; 
-V2(7 + 7r), n = l; 
L -V27, n > 2. 



Note that the it-solutions collide for 7 = tt. The linearization DF{^^ ^ 0) is invertible 
for 7 > TT, hence the implicit function theorem can bo applied again and we have 
the existence of monotone semi-kinks $o(^) + sf^^^kp-i^il)- In analogue with the 
continuum case, the type 1 wave is the one that has the discontinuity at the lowest 
value of the phase. The critical bias current for the existence of a static lattice 
semifluxon follows immediately from the arguments above. □ 

The two ib-solutions near 7cr as derived above in the proof are like the type 1 and 
type 3 semi-fluxons near 7cr in the PDEs studied in the previous two sections. So in 
analogue to those PDEs, we can define for 7 > tt 

$^(n;£;l-£7) = $o(n)+\/£$"(n;£;7) and $^(n;£;l-£7) = $o(«)+\/£5"^(n; £; 7)- 
where $0 and 6^ are as in the proof above. Thus we get 



(5.5) 



f - ^e ^ + 0{e ), n > -1; 

f TVi7 2(f^ + 0(£), n = 0; 

f - V^v/2(7 + ^) + 0(e), n=l; 

^/i^/2^ + +0{e), n > 2; 



2 

3tt 
2 

2 



The spectral stability of $^(n; e; 7) is obtained by substituting cpn = ^win; e; 7) -|- 
in the model equation (2.3). Disregarding the higher order terms in w„ gives 



At 



the following eigenvalue problem 



(5.6) 
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where A = A^, f = (. . . , V-i,Vo, Vi,. . .)'^ and L^{e; 7) is the linear discrete operator 

/•.•.•. n \ 



-2e-A_i e 

e —2e — Aq e 

e -2s - Ai 



A„ = cos (n; e; 7) + 6'„) , n e Z. 

This operator plays a similar role as the differential operator £'(2;; 7) = D^x — 
cos(05r(a;; 7) + in section 3. The eigenvalue problem is an infinite dimensional 

matrix problem for a real and symmetric matrix. Thus the eigenvalues must be real. 

In the discrete case, the continuous spectrum of semikinks is bounded. The 
spectrum is obtained by substituting w„ = g-tkn (5.6) with = —2e — — 7^ 
from which one obtains the following dispersion relation for such linear waves 

(5.7) A = -(v/r^ + 4£sin2(|)). 



Thus the continuous spectrum consists of the intervals ±i[\/l — 7^ , y \/l — 7^ + Ae] 
(recall that A = A^). 

In the following two lemmas we will show that all eigenvalues of the lineariza- 
tion L^{e;j) are negative for e small. Thus for e small, the type 1 wave is always 
stable. For 7 = 1 — £7 and 7 > tt, it can be shown that the linearization L^{e; 1 — e-y) 
has a positive eigenvalue for e small. Hence the type 3 wave is unstable for 7 near 7cr 
and with e small. 

Lemma 5.2. For < 7 < 1 and e small, the largest eigenvalue of the operator 
L^{£;"f) is negative up to 0{s'^). 

Proof. The eigenvalue problem to calculate the stability of the monotone discrete 
TT-kink <?i)^(n; e; 7), n e Z is given by (5.6) with i = 1. Slightly modifying Baesens, Kim, 
and MacKay [2] , the spatially decaying solution that corresponds to an eigenvalue of 
the above eigenvalue problem, can be approximated by 



(5-8) Vn = l^ 



c^-", n < 0, 
cc£"-\ n > 1, 



for some c, c and \£\ < 1. The diagonal elements in L^{e; 7) are An = — 7^+0(e^), 
if n 7^ 0, 1, Ao = - + 0{e^) and A^ = ^JY^ + e^=^ + 0{e''). 

Thus Aq ^ Ai, hence the need for two parameters c and c (modifying [2], where 
c = ±1 following from the symmetry Aq = A\). 

For small nonzero e, if we can match exponentially decaying solutions (5.8) on 
both sides from either end of the lattice to a central site, then we obtain a candidate 
for an eigenfunction. With (5.6), the parameters £ and c will be determined up to 
order e. For n ^ 0, 1, the relation (5.6) gives up to order e 



(5.9) A = -Vr^ + £(£-2+i). 
At the central sites n = 0, 1 we get up to order e 

(5.10) A=_^r3:^ + e^2^+e(i!-2 + c); 

(5.11) A=_^T^_,_2^+,(^_2 + i). 
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Combining (5.9), (5.10) and (5.11) shows that there are two possible value for c, being 

c± = /'^ ± V^+^^l^ — ^ )2_ _^ Oie) and leads to the eigenvalue A and the decay 

exponent £ as a function of e and 7, i.e. 



(5.13) K± = -^/l^+^{£±-lf + 0{e^). 

«± 

General Sturm-Liouville theory states that a critical eigenfunction that corre- 
sponds to the largest eigenvalue of a continuous eigenvalue problem does not vanish, 
except probably at x ^ ±00. This theorem can also be extended to a discrete eigen- 
value problem such that the most critical eigenvector does not have sign changes [1] . 
Thus, if we have a solution of the form (5.8) with £ > 0, then it is the critical eigen- 
vector. 

From (5.12), we see that 1+ > 0, thus the largest eigenvalue A+ from (5.13) is in 
the gap between zero and the interval associated with the continuous spectrum, i.e. 
A+ < 0. □ 

REMARK 5.3. From the details in the proof, note that weak coupling with strong 
bias current leads to one additional eigenvalue associated with ^- , where £_ < and 
1^-1 < 1. This indicates that the eigenvector of the form (5.8) is localized but has 
out of phase configuration, i.e. has infinitely many sign changes. This is a typical 
characteristic of a 'high-frequency' eigenvalue which is confirmed by the fact that A_ 
is indeed smaller than the phonon band. The presence of a high-frequency eigenvalue 
of a kink was previously reported by Braun, Kivshar, and Pcyrard [6] in their study 
on the Frenkel-Kontorova model with the Peyrard-Remoissenet potential [26] . 

For 7 close to 1, i.e 7 = 1 — ej, with 7 > tt, both the type 1 and the type 3 wave 
as given in (5.5) can be analyzed. In the following lemma we also show that both 
types have a high-frequency eigenvalue. 

Lemma 5.4. For 7=1 — ej, with 7 > tt, the largest eigenvalue of the opera- 
tor L^{e;l — ej) is strictly negative and the largest eigenvalue of the operator L^{e; 1 — 
£7) is strictly positive. 

Proof. As before, we write for an eigenfunction 

_ J c£-", n < 0, 
\ cc£"-\ n> 1, 

for some c, c and \£\ < 1 and we substitute this in the eigenvalue problem, which leads 
to the equations 

(5.14) A = - sin + C'(e)) +e{l/i-2+ 1)- 

(5.15) A = T sin (^2^(7 - 7^) + 0(e)) + e(f - 2 + c); 

(5.16) A = - sin (V2£(7 + 7r) + 0{e)) + e(f - 2 + 1/c). 

where =p-sign in the second equation is a minus sign for the eigenvalue problem asso- 
ciated with the type 1 wave and the other sign in case of the type 3 wave. Again, by 
subtracting (5.16) from (5.15), we get a quadratic equation for c, with two solutions. 
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one of order and one of order yje (which is easiest found by writing the equation 
as a quadratic equation for 4): 



j_ _ 



C2 



7? (- v/2(7 + tt) ± ^2(7 - tt) + O(V^) . 



Combining (5.14) and (5.15) rcsp. (5.16) and using the two expressions above gives 
that in both cases I is of order \/e and given by 



I = ^fyWT\/2(7-^) + 0(V^) 



Finally, substitution into (5.14) shows that 



(5.17) 



Ai = T\/g V2(7-7r) + 0(£); 
A2 = -^/iV2(7 + 7^) + 0(£). 



The eigenvalue that corresponds to £ > is Ai. 

So clearly the largest eigenvalue Ai is negative in case of the type 1 wave and is 
positive in case of the type 3 wave. 

In addition, the operator L^{e; 1 — £7) and L^{e;l — £7) have the same high- 
frequency eigenvalue A2 (up to order e^). D 

The proofs of lemmas 5.2 and 5.4 show the presence of a high-frequency eigenvalue 
for a semi-kink in case the bias current is not small. In the following we will show 
that the eigenvalue appears when the bias current is larger than + 0{e). 

Because we do not have an analytic expression for the type 2 and type 3 semi-kink 
in the small forcing limit, the analysis is done only for the type 1 semi-kink. 

Lemma 5.5. There is a critical value 7hf, 7hf = \/e + 0(e) such that for all 
7 € (7hf,7cr) the operator L^{e;j) has a high frequency eigenvalue that up to 0{e^) 
is attached to the Imnest boundary of the continuous spectrum. The corresponding 
eigenvector is localized and changes sign between any two adjacent sites. 

The appearance of this eigenvalue and the structure of its eigenvector is checked 
numerically in section 6. 

Proof. Again, we write for an eigenvector 

c^-", n < 0, 
cc£"-\ n > 1, 

for some c, c and \C.\ < 1 and we substitute this in the eigenvalue problem. We 
first consider 7 = ^£7- This gives An = l-^-^+ Oie"^'^), if n ^ 0, 1, = 
l-^-£3/M-£'(^ + ^)+^?(e'/')andAi = l-^+£3/M-£'(^ +^)+g(g'/')- 

Using the same procedures, this implies c± = —^/en^—e^/'^'y^n±^/^^^^^F+'2^T^^^e^'^\^ 

and 

^ = £77^ ± 2^7727% + 27727^^2 + 1 = ±1 + £^"(^^7') + 0(g2-)_ 

l± 2 

For 7 > 1 there are two solutions \£±\ < 1; > corresponds to the largest 
eigenvalue and also exists for 7 < 1 (Lemma 5.2); £- < 0, so that its associated 
eigenvector w„ indeed changes sign between any two adjacent sites. 
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Fig. 6.1. (a) Two lattice sernifluxons of type 1 with no bias current (7 = Oj are plotted as 
a function of the lattice index, namely the kink for strong coupling with e = 100 (equivalently, a 
very small lattice spacing a = O.l) (— * i.e. close to (3.10) and tlie kink for weak coupling with 
e = ^ (equivalently, a large lattice spacing a = 2) (—o—). (b) Numerically computed spectrum of a 
lattice semifiuxon against the lattice spacing parameter a with 7 = 0. We used the number of sites 
2N = 300. We zoom in the plot of spectra around -1 for clarity. The bold-solid-line is the calculated 
approximate function for the point spectrum using perturbation theory for a small resp. for e small. 

The value 7hf = ^/s-\-0{e) indicates the appearance of this high frequency eigen- 
value from the continuous spectrum. It follows from a straightforward analysis that 
this eigenvalue exists, i.e. £- € (—1,0) exists, for all 7 G (7hf,7cr), and the corre- 
sponding eigenvalue is A_ = -1 + 4)e + ^"^/Se^ + 0{e^). Up to 0{e^) this 

eigenvalue is nothing else but the lower boundary of the continuous spectrum. □ 

6. Numerical computations of the discrete system. To accompany our an- 
alytical results, we have used numerical calculations. For that purpose, we have made 
a continuation program based on Newton iteration technique to obtain the stationary 
kink equilibria of (2.3) and (2.4) and an eigenvalue problem solver in MATLAB. To 
start the iteration, one can choose either the continuum solutions discussed in sec- 
tion 3, i.e., the case where the lattice spacing parameter a = 0, or trace the equilibria 
from the uncoupled limit e = (a ^ 00) as discussed in the previous section. We 
use the number of computational sites 2N = 800 for parameter values of a = 0.05 or 
larger (e = 20 or lower). 

6.1. Stability of type 1 lattice semifiuxon. The type 1 lattice semifiuxon 

<I>^(n;£;0), n G Z has been studied analytically both in the strong coupling limit 
(a <C 1, or e » 1) and the weak coupling limit {a ^ 1, e 1). In Figure 6.1(a), 
$^(n; e; 0) is plotted for two different values of the coupling parameter e. For a given 
value of e, one can use as initial guess in the numerical procedure either a solution 
from the continuous limit (3.10) or from the uncoupled limit that has been discussed 
in the preceding sections. 

In Figure 6.1(b), we present the numerically calculated spectrum of the type 1 
semifiuxon with 7 = as a function of the lattice spacing parameter. The approximate 
largest eigenvalue (4.9), derived for a small, and the one derived in Lemma 5.2 for a 
large, are in a good agreement with the numerically obtained largest eigenvalue. Any 
eigenvalue below A = — 1 belongs to the continuous spectrum. For a close to zero we 
do not see dense spectra because of the number of sites we used. By increasing the 
sites-number we will obtain a more dense spectrum. 

There is only one eigenvalue outside the phonon bands, the largest eigenvalue 
as studied in Lemma 5.2. This is in contrast to the case of an ordinary lattice 27r- 
kink [17, 19] where there is an internal mode bifurcating from the phonon band when 
the parameter a increases. 

If Fig. 6.1(b) shows the spectrum of the type 1 semifiuxon as a function of the 
coupling parameter e {e = 1/a^) for a fixed bias current 7, in Fig. 6.2 we present the 
numerically calculated spectrum of the type 1 lattice semifiuxon as a function of 7 for 
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Fk;. 6.2. (a) Spectrum of the type 1 semifluxon as a function of the applied bias current 7 for 
a value of the coupling constant e = 0.25. The dashed line is theoretical prediction from (5.13). 
In (b) we zoom in on the spectrum around —1 for clarity. The spectrum is normalized to the lower 
edge of the phonon band, i.e. ^1 — 7^ +4£, such that the appearance of a high frequency eigenvalue 
can be seen clearly. 



Fig. 6.3. The eigenvectors associated the two largest eigenvalues and the two smallest eigen- 
values of the truncated 2N x 2N -matrix associated with L^{e;'y) (the type 1 discrete semikink) for 
e = 0.25 and 7 = 0. The results are shown with 2N = 100 for clarity. Shown are the eigenvector 
of (a) the largest eigenvalue, (b) the second one, (c)-(d) the last two eigenvalues. There is only one 
eigenvalue for (e; 7) as there is only one localized eigenvector. 



Fig. 6.4. The same as Fig. 6.3 for 'y = 0.5. 



a fixed e, e = 0.25. 

Lemmas 5.2 and 5.5 established the existence of two eigenvalues (for e small 
enough) for the stability problem associated to the type 1 semifluxon, the largest 
eigenvalue and an additional eigenvalue which bifurcates from the lower edge of the 
phonon band for bias current 7 > 7ht. It follows from the numerical simulations that 
these arc indeed the only two eigenvalues (Figure 6.2). For e = 0.25, this inininnun 
bias current 7hf is approximately 0.466. Interestingly, according to Lemma 5.5 the 
bifurcation appears at 7hf = \/e = 0.5 at leading order in e. This is in remarkably good 
agreement with the numerical result, especially since the error is 0{e) and e = 0.25. 

To picture the appearance of the high-frequency eigenvalue, all eigenvalues for the 
truncated 2N x 2iV-matrix associated with L^(e;7) are determined and the eigenvec- 
tors of the two largest eigenvalues and the two smallest eigenvalues are presented in 
Figures 6.3-6.5 for various values of 7 and a fixed e. It can be observed that there is 
always a localized eigenfunction associated with the largest eigenvalue. In Figure 6.3 
(7 = 0), none of the other eigenvectors can be associated with localized eigenfunctions 
and in Figures 6.4 and 6.5, the birth of the localized eigenfunction associated with 
the smallest eigenvalue can be observed. 

If we keep increasing 7 further, then there is a critical applied bias current at 
which the largest eigenvalue becomes 0. Numerical computations show that this 
critical value is 7cr(£) above which static lattice semifluxons disappear. 
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Fig. 6.5. The same as Figure 6.3 for ^' = 0.7. Note that there are now two localized eigenvectors 
shown in (a) and (d). The smallest eigenvalue associated with (d) is —1.7357 while the lower edge 
of the phonon band is —1.7141. Note also that neighboring sites of the eigenvector in (d) move out 
of phase indicating a high-frequency mode, contrary to the semikink's low-frequency mode in (a). 



Fig. 6.6. The critical bias current of a static ir-kink as a function of the lattice spacing pa- 
rameter a. For 7 above the critical current there is no static n-kink solution. The solid line is 
numerically obtained curve. Dashed lines are the theoretical predictions (4-^) for o <C 1 resp. (5.4) 
fora:S>l (e<^ 1). 



Fig. 6.7. Plot of the eigenvalues of a Sir-kink as a function of tlie lattice spacing parameter a. 
We zoom in tfie rr.gion with a <g; 1 where it shows that turning the lattice spacing on destabilizes 
the kink. The dashed line depicts the analytically computed approximation (4-10) to the largest 
eigenvalue of the 3n-kink. 



The critical bias current for the existence of a static type 1 lattice scmifluxon in the 
continuum limit and for a very weak coupling in the discrete system has been discussed 
and analytical expressions were given in sections 3, 4 and 5, respectively. In Figure 6.6, 
the numerically calculated critical bias current 7cr of the discrete system (2.3) as a 
function of the lattice spacing a is presented. The approximate functions, given 
in (4.2) for small a, and in (5.4) for large a (small s) are presented as dashed lines. 



6.2. Instability of type 2 lattice semifluxon. In the continuum models we 
have seen that for 7 small, the instability of the type 2 semikink is mainly determined 
by the instability of the 37r-kink in the continuum models for 7 = 0. So we start 
this section by looking at the stability of the 37r-kink in the discrete model. We will 
denote the 37r-kink by ^"^^^{11; s; 0), where as before the coupling parameter e and the 
lattice spacing a are related by e = ^ . 

Using our continuation program, we have followed a 37r-kink solution from the 
continuous limit < o <C 1 up to the uncoupled situation £ = (i.e., a = oo). We 
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Fig. 6.8. (a) Plot of a type 2 sernikink with 7 = 0.01 for e = 100 (- * -) and e = 40 (-0-). 
(b) Plot of the largest eigenvalue of a type 2 sernikink as a function of the coupling parameter e. 
When e = 0, the eigenvalue converges to A = \J\ — . 



obtain that '^\^(n\ 0; 0) is given by 



(6.1) ^lAn\m 



0, n = -l,-2, . 
27r, n = 0, 
TT, n = 1, 
37r, n = 2,3, ... 



Note that this discrete configuration is not monotonically increasing as opposed to 

the continuum configuration, which is monotonic. 

In Figure 6.7, we present the numerically obtained eigenvalues of a 37r-kink as a 
function of the lattice spacing a. For small a, the largest eigenvalue is indeed increasing 

as is predicted by the perturbation theory (4.10). As soon as the lattice spacing is 
of order one, the largest eigenvalue decreases and becomes zero at approximately 
a = 1.7521. 

After establishing that increasing the lattice spacing can stabilize the 37r-kink at 
7 = 0, we continue by looking at the stability of the 27r-kink for 7 > 0. Interestingly, 
increasing the lattice spacing does not stabilize a type 2 semikink for 7 > 0. In 
Figure 6.8, we show a plot of the type 2 semikinks for two values of e as well as a 
plot of the largest eigenvalue as a function of e for two particular values of 7, namely 
7 = 0.01 and 7 = 0.1. We present the largest eigenvalue as a function of the coupling e 
instead of the lattice spacing a as the eigenvalue changes most for small coupling (large 
lattice spacing). From Figure 6.8 it follows that the solutions are unstable even in the 
weak-coupling limit. This is interesting as in the limit for 7^0, the type 2 semikink 
can be seen as a concatenation of a 37r-kink and a — 27r-kink. Both the 37r-kink and 
— 27r-kink arc stable for the coupling e sufficiently small, while the type 2 semikink 
turns out to be unstable. 

This instability issue can be explained by looking at the expression of a type 2 
semikink when it is uncoupled (e = 0). For the two particular choices of 7 above, we 
get from the simulations that the configurations of the these semi-kinks are given by 



(6.2) ^l{n;^-^m) = 



+ arcsin(O.Ol), n < -1, 

TT — arcsin(O.Ol), n = 0, 

TT + arcsin(O.Ol), n = l, 

37r + arcsin(O.Ol), 2 < n > 8, 

27r — arcsin(O.Ol), n = 9, 

TT -|- arcsin(O.Ol), n > 10, 
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Fig. 6.9. (a) Plot of a type 3 semikink with 7 = 0.01 for e = 100 (- * -) and £ = 40 (-0-). 
(b) Plot of the largest eigenvalue of a type 3 semikink as a function of the coupling parameter e. 
When e = 0, the eigenvalue converges to A = \J\ — . 



and 



(6.3) 



$2(n;0;0.1) 



+ arcsin(O.l), 
TT — arcsin(O.l), 
TT + arcsin(O.l), 
Stt + arcsin(O.l), 
27r — arcsin(O.l), 
TT + arcsin(O.l), 



n < -1, 
n = 0, 
n = 1, 
2 < n > 6, 
n = 7, 
n > 8. 



We see that there are two sites, namely n = and n = 9 for 7 = 0.01 and n = 
and n = 7 for 7 = 0.1, where <I> takes the vahic of an unstable fixed point of the 
discrete system (2.3). Looking only at sites numbered n = 2 to n — + 00, $^(n;0;7) 
can be viewed as a — 27r lattice kink sitting on a site which is known to be unstable. 
If wo look only at sites numbered n = 6 to n ^ —00, <l>^(n;0;7) can be seen as a 
deformed 37r lattice kink (at site n = 0, the phase $ takes the value n instead of the 
value 27r as in the 37r-kink). Hence, it seems that coupling between the two kinks due 
to the presence of a nonzero 7 is responsible for the instability. 

It has been discussed in the previous sections that there is a critical bias current 
7* for the existence of a type 2 lattice semikink in the continuum models. However, we 
did not numerically calculate the critical bias current 7* (a) for discrete system (2.3). 



Fig. 6.10. Spectrum of the type 3 semifluxon as a function of the applied bias current 7 for a 
value of the coupling constant £ = 0.25 . In (b) and (c) we zoom in near the phonon band for clarity. 
In (b), the spectrum is normalized to the upper edge of the phonon band, i.e. -y/l — 7^ and in (c) 
it is normalized to the lower edge of the phonon band, i.e. \/l — 7^ + 4e. The disappearance of a 
high-frequency mode in the lower edge of the phonon band can be clearly observed in (c). The insets 
in (b) and (c) show the eigenfunctions of the two eigenvalues just above resp. just below the phonon 
band for 7 = 0.73. 

6.3. Instability of type 3 lattice semifiuxon. In this section, wc will consider 
the type 3 semikinks, which will be denoted by $^(n; s\ 7). In Lemmas 3.9 and 4.7 it 
has been shown that these kinks are unstable in the continuum models for small or 

zero lattice spacing. 

The largest eigenvalue of a lattice type 3 semifluxon for three particular values 
of 7, i.e. 7 = 0.01, 0.1, 0.55, is presented in Figure 6.9. Even though in the limit 
for 7 — > 0, a semifluxon of this type is a concatenation of a 27r-kink and a — 7r-kink 
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which both can be stable in the discrete system, the type 3 semikink is unstable for all 
parameter values from the zero lattice spacing limit all the way to the zero coupling 
one. The explanation is similar to the one for a type 2 semikink discussed above. 
Indeed, for the three particular choices of 7 above, 0; 7) is given by 

n = -l,-2,... 
n = —6, 
n = —5, . . . , 0, 
n = 1,2, 



= -1,-2,... 
= -2, 
= -1,0, 
= 1,2,..., 



= -1,-2,... 

= 0, 

= 1,2,.... 

One interesting point to note for the type 3 semikink is that the number of sites with 
value 2tt is decreasing as 7 increases. Starting from the continuum approximation of 
a type 3 semikink as the initial guess for the continuation program, the 27r-plateau 
disappears for 7 > 7*(a) (see (4.1)). For 7 > 7* the configuration at e = is similar 
to the stable typo 1 7r-kink (5.2), apart from the value of the phase at the site with 
n = (where the phase takes the value of an unstable fixed point). 

Because analytical calculation of the spectrum of the type 3 semikink has been 
obtained in the small coupling Umit and bias current close to 1 (Eqs. (5.17)), it is worth 
comparing the analytical predictions with numerical computations. The theoretical 
calculations shows that for e small and 7 close to 1, the type 3 semikink has at least 
two eigenvalues, one of which corresponds to a high frequency mode and the other to 
a positive eigenvalue. 

Using the continuation of (6.6) for e = 0.25, the spectrum of the type 3 lattice 
semikink is presented in Figure 6.10 as a function of the applied bias current 7. 
Our numerics show that when 7 is very close to 7cr, the type 3 semikink has three 
eigenvalues, one of which corresponds to a high frequency mode and is below the 
phonon band, while the other two are above the phonon band. The birth of this high- 
frequency mode is shown in Figure 6.10(c) and is qualitatively similar to the case of 
the type 1 lattice semikink. The two eigenvalues which exist for all values of 7 can 
be observed in Figure 6.10(a) and (b). 

7. Conclusions. We have performed an existence and stability analysis for three 
types of lattice 7r-kink solutions of the discrete O-vr sine-Gordon equation and its 
continuum limits. Analytical results have been established in the continuum limits 
and in the weak-coupling case. It has been shown that in the continuous O-tt sine- 
Gordon equation, 7r-kinks of type 1 are stable and the other types are unstable. 
The introduction of discreteness destabilizes the unstable 7r-kinks even more. An 
approximation to the largest eigenvalue of all types of 7r-kinks has been derived both 
in the continuum and the weak coupling limits. 



0-f arcsin(O.Ol), 

,\ ^ I TT — arcsin(O.Ol), 

(6.4) $^(n;0;0.01) = { r, , • /n m\ 

^ ' TTv I I y I 27r + arcsm(O.Ol), 

T -|- arcsin(O.Ol), 



-|- arcsin(O.l), n 

rx « I TT — arcsin(O.l), n 

(6.5) $^(n;0;0.1 = { „ , - Vn i\ 
^ ' 7r\ ' ' / j 27r -I- arcsm(O.l), n 

TT -|- arcsin(O.l), n 
and 

-I- arcsin(0.55), ? 

(6.6) ^l{n;0;0.55) = { tt - arcsin(0.55), r 
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For future research, it is of interest to study the nucleation of kinks and antikinks 
when a constant force, or bias current, 7 that is above the critical value j^r is appUed 
- sec Figure 3.2(b). One question that can be addressed is the mechanism and the 
frequency of the nucleation as a function of the applied constant force, especially in 
the presence of a damping coefficient (which has not been considered in this paper). 
In work in progress, the stability of the type 3 semifluxons in the presence of defects 
in studied. These semifluxons are unstable, but the largest eigenvalue is close to zero. 
In fact, a type 3 scniifluxon consists of a fluxon and a semifluxon with the opposite 
polarity. In experiments, the presence of a fluxon nearby a semifluxon can influence 
a junction measurement [12]. Because a fluxon can be pinned by a defect [18], one 
can expect to have a stable type 3 semifluxon when there is a defect present in the 
system. 

Acknowledgments. H.S. wishes to thank Panayotis Kevrekidis for numerous 
useful interactions and discussions. 
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